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ABSTRACT 


A parametric study of pressurized light water reactor 
cores was undertaken, in which both heat transfer and nuclear 
influences of the basic parameters on the quantity and required 
enrichment of the uranium fuel were to be considered. 


Methods were developed for the following: 

Ll. Computation of the minimum core size for a specified 
power output on the basis of heat transfer characteristics of 
the unit cell parameters. 

2- Computation of the critical mass of U-235 required 
for a given reactor by an adaptation of a digital computer 
program written by J. R. Powell. This program, written for 
Whirlwind, is a Fourier solution for critical mass of a highly 
enriched, homogeneous, heavy water moderated reactor; methods 
were developed to adapt it to the low enrichment, heterogeneous, 
light water case, and to compute the necessary inputs. 

3. Computation of the factor by which the critical en- 
richment of a given reactor would have to be increased in order 
to attain a specified core life. 


A correlation was obtained between the adjusted Powell 
Fourier program and experimental data; the methods selected 
for calculation of resonance escape probability, fast fission 
factor, and disadvantage factor for use in the program were 
also confirmed by comparison of results with experimental data. 


Ranges of parametric variation were chosen as follows: 
lL. Fuel elements to be rods, 0.2, 0.4, and 0.6 inches 
in diameter. 
Ze Ratio of volume of moderator to volume of fuel to 
me 1.5, 2.0,°3.0, 4.0, 50, dad 6.0. 
3- Cladding to be of stainless steel and of zirconiun, 
0.015 inches thick. 
4, Fuel to be uranium oxide, UOo; quantity and enrich- 
ment to be determined. 
5. Cores to be right circular cylinders not to exceed 
13 feet in diameter and core heights to be 4, 6, and 8 feet. 
6. Cores to be surrounded by an effectively infinite 
water reflector. 
7. Maximum fuel temperature: 4800°F. - 
Coolant pressure: 2000 psia. Inlet temperature: 487°F 
Maximum t emperature: 62 OF 
Average reactor temperature: 5080°F 
8, Thermal power output: 400 megawatts. 


For each case, an appropriate core diameter and total 
quantity of uranium were computed from the heat transfer charac- 
teristics associated with the parameters. With this, resoriance 
escape probability, fast fission factor, and cross sections 
modified for disadvantage factor, the machine calculation was 
entered. The output of the machine calculation was critical 
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mass of U-235; this was compared with the total quantity of 
uranium previously computed to give the required critical en- 
richment. 


Analysis of the uranium quantity and enrichment data 
thus obtained indicated that: 

lL. Total quantity of uranium required, a function of 
the heat transfer characteristics of the unit cell, is ap- 
proximately proportional to the square of the rod diameter; 
the allowable heat generation rate per unit length of the cen- 
tral rod is very nearly the same regardless of the rod diam— 
eter, and hence the number of rods required is very nearly 
the same for all rod diameters at a given height. 

2. Critical enrichment is principally a function of the 
unit cell parameters and is nearly independent of the core 
hae and shape between diameter to height of core ratios of 
0.8 and 3. 

3. For each rod diameter and cladding material an opti- 
mum moderator to fuel volume ratio exists at which the criti- 
cal enrichment is a minimum. For the diameters tested, these 
were as follows: 


Optimum V1/y> 


Rod 
ae Dione vem Stainless Steel clad Zirconium clad 


0.6 inch 3.9 3,0 
0.4 inch 3.2.75 4.0 
0.2 inch 45 5 


4, An optimum rod diameter at which minimum critical. 
enrichment is lowest exists for rods clad with each eer ga 
material. For zirconium clad rods, this occurs at about 0.45 
inches rod diameter for stainless steel clad, at slightly 
over 0.6 inches. 


Rod Diameter Minimum critical enrichment 
Stainless Steel clad Zirconium clad 

0.6 ineh 99 89 

0.4 inch 1.12 ~85 

0.2 inch 1.5 ae 15 


5. Critical enrichments are lower for zirconium clad 
rods than for stainless steel clad, of the order of 1% en- 
richment compared to the order of 24. 

6. Optimum critical enrichment does not necessarily in- | 
dicate the optimum core; other factors which must be considered 
are total quantity of uranium, cost data for uranium, cladding 
materials, and fabrication of fuel elements and pressure ves-— 
sels, and the additional enrichment necessary to attain a 
specified core life. 
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NOMENCLATURE 


A. General 


Symbol 


Description Dimensions 
Atomic weight 
Coolant channel cross sectional area ft 
Rod heat transfer area ft 
Buckling-a measure of the bending em 7 
of the neutron flux at a point in 
a critical reactor 
Coolant specific heat BTU/Lb-F 
Diffusion coefficient-the proportion- 
ality factor between neutron current 
and neutron flux gradient 
Energy Mev 
Dimensionless coefficient=-related to 
the computation of the disadvantage 


factor. <A function of Kand R. See 
Appendix A. 


Hot channel factor-heat transfer 
safety factor 


Thermal utilization-the proportion of 

thermal neutrons absorbed in fuel to 

the total absorbed 

Core height cm 

Film heat transfer coefficient BTU/hr-ft2-P 


Modified Bessel function of the first 
kind 


Bessel function of the first kind 


Modified Bessel function of the second 
kind 









Infinite multiplication factor-the 
ratio of the average number of neutrons 
produced in each generation to the 


average number of corresponding neu- 
trons absorbed 


Thermal conductivity 


(Diffusion length)*=-1/6th of the mean 


square distance traveled by a thermal 
neutron before being absorbed 


Migration area-1/6th of the mean 
square distance from birth of fission 
neutron to capture of a thermal neutron 


6.025 x 10°2 
—_ noel 


Dimensionless coefficient-function 
related to temperature differences in 
heat transfer. See Appendix G. 


Resonance escape probability- the 
probability that a neutron will not be 


absorbed by U-238 in the resonance 
energy region 


Dimensionless coefficient-function 
related to heat transfer resistance 


Slowing down density-rate at which 
neutrons reach thermal energy per 
unit volume 

Heat transfer rate. See Appendix G. 
Radius 

Temperature 


Total heat transfer coefficient 


Volume 


Ratio of volume of uranium 4 clad to 
water in a unit cell 


Fast fission factor-the proportion of 
total fissions at both fast and thermal 
energies in both U=235 and U=-238 to the 
number of thermal fissions in U-235 


BTU/hr-ft-or 


om 


Ee 
cm 


atoms/ em? 


neutrons/ 
cm--sec 


BTU/hr-ft~ 


Op 
BTU/hr-ft2.Or 


em3 


iv 





= 2 


Etr 


Reciprocal diffusion length-~ = 1/L 


Reflector savings-the decrease in 
critical size of a reactor due to 


a reflector 


Density 


Microscopic absorption cross section- 
the probability that one atom of 


material will absorb a neutron 


Microscopic fission cross section- 
the probability that one atom of 
material will be fissioned by a 
neutron 


Microscopic scattering cross section- 
the probability that one atom of 
material will scatter a neutron 
Microscopic transport cross section-=- 
the product of ce and (1- Ao) 


Macroscopic absorption cross section- 
the product of the number of atoms of 


material per unit volume and G”’for tha 


material | 


Macroscopic fission cross section= 


the product of the number of atoms 
of material per unit volume and 7 
for that material 


Macroscopic transport cross section- 
the product of the number of atoms of 


material per unit volume andG@y for 
that material 


Age-1/6th of the mean square distance 
traveled by a neutron in being slowed 
down from fission to thermal energy 


Average number of fast neutrons re- 
leased per slow neutron fission 


Neutron flux-the product of the 
number of neutrons per unit volume 
and their velocity 


barns 


barns 


barns 


em~1 


em~L 


em~1L 


ome 


neutrons/ 
em’-sec 
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B. Specific 


Symbol 


9 


ag 


Description 


rod heat transfer area, midpoint to end of core 


rod heat transfer area midpoint extended through 
reflector savings 


eross sectional area of unit cell 

heat generation rate hot channel factor 
temperature rise hot channel factor 

heat transfer coefficient hot channel factor 
clad thermal conductivity 

fuel thermal conductivity 


reactor thermal power 


average power developed per rod averaged over core 


average power developed in a rod averaged longitudi- 
nally only 


heat transfer rate at midpoint of the central rod 


heat transfer rate at midpoint of the central rod 
required to produce burnout of clad 


heat transfer rate at midpoint of the central rod 
required to produce local boiling at outlet 


radius of fuel in a unit cell or to inside of cladding 
outside radius of unit cell 

radius to outside of cladding 

core radius 


core radius plus reflector savings 
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Ate 
A tr 
Atwi 


A two 


temperature rise of coolant in core from inlet to 
maximum temperature 


temperature rise from coolant inlet to centerline 
of fuel rod 


temperature rise from coolant inlet to cladding- 
fuel interface 


temperature rise from coolant inlet to outside of 
cladding 


effective heat transfer coefficient for bond, clad, 


and scale 
volume of fuel in unit cell 
volume of clad + moderator in unit cell 


volume of clad in wnit cell 


volume of water in unit cell 

volume of uranium in unit cell 

volume of structure and clad in unit cell 
fuel reciprocal diffusion length 

moderator reciprocal diffusion length 

fuel macroscopic absorption cross section 
moderator macroscopic absorption cross section 
fuel macroscopic transport cross section 
moderator macroscopic transport cross section 
slowing down cross section 

Fermi age 

diffusion age 

volume average slow neutron flux in fel 
volume average slow neutron flux in moderator 
volume average slow neutron flux in clad 


disadvantage factor 


vil 





oe 
INTRODUCTION 





The category of nuclear reactors known as pressurized 
water reactors has emerged from the multitude of possible types 
of power producing reactors as one of the most promising from 
a feasibility viewpoint. It appears that over the next sev- 
eral years it will continue to be an important type both for 
ship propulsion and for shore power station use, and hence a 
study of the effects of varying some of the basic core para- 
meters was considered to be a fruitful undertaking. 

A pressurized water reactor is distinguished by the fol- 
lowing characteristics: 

1. Water is used as both the coolant and the moderator. 

2. It is heterogeneous; that is, it is made up of 
segregated regions of fuel and coolant-moderator, 
separated by a corrosion resistant "cladding." 

3. The water is kept under high pressure to prevent 
boiling. 

4, The great majority of fissions is produced by 
neutrons at thermal energy, that is, neutrons 
which are in thermal equilibrium with the atoms 
of the surrounding medium. 

The first full seale power reactor to operate, the pro- 
totype of the propulsion reactor in the U. 8S. submarine 
Nautilus, was of this type. Since then, in addition to the 





successful plant actually installed in the Nautilus, the 
pressurized water concept has been accepted as standard for 


U. 8. submarines and thirteen other boats are building, in 
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the design state, or authorized. Aircraft carrier and cruis- 
er plants are also being planned. The pressurized water 
reactor has also taken a leading position in the U. 8S. for 
shore power station use: the first civilian power reactor 

to operate will be the PWR at Shipping port, Pennsylvania, 

and among the ones to follow it will be the Yankee Atomic 
Electric Company's design, to be built at Rowe, Massachusetts. 
Similar reactors are being built in other countries; a list 


hs. 
of those now planned follows:~— 


Heat 
Fuel Power Critical Date 
Army Package Power Reactor U05 10 mw Early 1957 
Belgian Thermal Reactor U 260 1957 
Dominican Republic Reactor 05 43 
Consolidated Edison Reactor U 500 July 1959 
University of Florida Reactor U 10 December 1959 


The uranium fuel in all these designs is formed into 
elements, such as plates or rods, which are small compared to 
the size of the reactor as a whole but very large compared to 
atomic distances and of the same order of magnitude as the 
mean free path of a neutron in uranium, about half an inch. 
These elements are regularly spaced in a vessel through which 
water can flow, filling the spaces between the elements. Each 
element with Lee mnaes of the surrounding water can be thought 
of as a cell, and the whole reactor as an aggregate of such 
cells, or a cell lattice. Some of the nuclear and thermal 


properties of the reactor are properties of the cell parameters 


1."Nuclear Reactor Data Handbook No. 2", The Raytheon Mfg. Co. 





while some are properties of the lattice. The properties of 


basic interest are —— 


l. Nuclear properties of the unit cell. 

a. The infinite multiplication factor ko. In 
order to maintain steady state operation, the 
number of neutrons in successive generations 
must be conserved. That is: 

Production = Absorption t Leakage 
Koo is the number of thermal neutrons produced 
per thermal neutron absorbed if no leakage 
takes place. It can be thought of as the 
factor relating the numbers of neutrons in 
successive generations in a medium consisting 
of an infinite repetition of the unit cell. 
It is generally recognized as the product of 
four factors, which, with the energy cycle 


followed by the neutrons, are explained below. 


Ul > The number of fast neutrons produced 
per thermal neutron absorbed in the 
fuel. High speed neutrons with a mean 
energy of about 1 mev are produced by 
absorption of thermal neutrons, those 


which have been so slowed down by 


L. These properties are described in more detail in various 
references, and, in some cases, in the Appendices to this 
thesis. These references are cited in each case. 

2. See §. Glasstone and M.C. Edlund, “The Elements of Nuclear 
Reactor Theory", Van Nostrand (1952)(hereafter cited as 
"Glasstone and Edlund"), Section 4.57 ff. 





successive collisons that their energies 
depend only on the temperature of the 
surrounding medium. Not all those absorbed 
in the fuel produce fissions; those that 
do, produce on the average 2.46 fast neu- 


trons and about 192 mev of energy. 


e : The fast fission factor. Neutrons with 
energies greater than about 1 mev are 
capable of producing a small number of 
fissions in U-238. g¢€ is the ratio of the 
total rate of neutron production to the 
production rate from thermal fission of 
U=235. 

The resonance escape probability. U-238 
has a high probability of absorption with- 


out fission for neutrons in the "resonance" 
energy range, from about 10 to 100 ev. p 
is the probability that a fission neutron 
in slowing down, will not be so captured. 

f : The thermal utilization is the proportion 


of all the thermal neutrons absorbed 
which are absorbed in the fuel. 


b. The disadvantage factor, += “fy, - Thermal 


neutrons are produced in the moderator region 


as fast neutrons resulting from fissions in 


ee 


1. See Glasstone and Edlund, 9.48, and Appendix A. Also see 
A. L. Kaplan, “Theoretical Studies of Neutron Flux Distri- 
bution Expected in the MIT Nuclear Research Reactor", MIT 
S.M. thesis, 1955. 





the fuel region are slowed down by elastic 
collisions with the hydrogen nuclei. 

The neutron flux = p , expresses the length 
of neutron paths covered per unit time per unit 
volume. This depends both on the number of 
neutrons and on their velocity. This is shown 


in the units of flux: 


neutrons x — x = = age 
; cm om~“=gec 


The number of interactions taking place per 
unit time and volume is the product of the 
length of neutron paths and the probability 


of interaction per unit path length, S,. 


no. of interactions - Ss p 
unit time=unit volume 


= in turn consists of the number of nuclei 


per unit volume, N, times the reaction proba- 


bility per nucleus,o”. Thus, 
no. of interactions = Zo = No- 


unit time~unit volume 

The fast flux produced by fissions in the 
fuel is slowed down in the moderator, producing 
thermal fiux. As the thermal flux thus produced 
reenters the fuel, it is continuously reduced 
as the center of the fuel element is approached 
since & for absorption is greater in the fuel 
than in the moderator and since essentially 


no slowing dow takes place in the fuel to 


1. See Glasstone and Edlund, 3.49. 





replenish the supply of thermal neutrons. The 
general pattern of thermal flux in a cell is 
shown in Figure l. 

The fuel near the center of the element, 
thus, meets a lower flux than the average over 
the whole reactor, with a consequent reduction 
of reaction probability below that which would 
be exhibited by a homogeneous mixture of the 
same materials in the same average thermal 
flux. The ratio of the average thermal neutron 
flux in the moderator to that in the fuel is 


the disadvantage factor. 





Ne, 
. Fermi Age, Q& . Since the slowing down process 


is a complex one, involving fission neutrons 
with a distribution of initial energies taking 
part in a large number of elastic collisions 

in which their energy loss varies (for a hydro- 
gen containing medium) between zero and 100%, 
it is represented by an experimentally deter- 


mined number which is characteristic of the 


1. See Glasstone and Edlund, 6.117 to 6.145, and Appendix F. 





proportions of the various materials in the 
reactor; since these materials and proportions 
are repeated in all the cells, it can thus be 
considered a property of the cell. Physically, 
the Fermi age is one sixth of the average 
squared crow-flight distance travelled by a 
neutron during the process of slowing down 
from fission to thermal energy. 

2. Nuclear properties of the reactor as a whole. 

ae Flux distribution, = P+). Due to the effects 

of leakage and changes of composition across 
the lattice, the flux varies over each dimen- 
sion of the reactor, being high near the center 
and in regions of high multiplication factor, 
and low near the edges and in regions of low 
multiplication factor. It drops to zero a 
small distance outside the physical boundaries 
of the reactor at the “extrapolated poutiaawy, *=— 
The functions which express these patterns are 


as follows: 


Reactor Shape Dimension Function 
Cube parallel to an edge ear x 
of length H 2 4/2 


Right cirecu- radial(radius = R) Jo (2-505n) 
lar cylinder R 


axial (height = H) c08(E Fe) 


Sphere radial (radius = R) Te sin(Tr ) 
R 


I. See Glasstone and Edlund, 7.36 to 7.61 
2. See Glasstone and Edlund, 5.40 














A reflector is frequently provided to reflect 
neutrons back into the core which would other- 
wise leak from it. When this is done, the 
average flux in the active core is increased, 
as shown in Figure 2, and the flux pattern 

is extended to the reflector boundary, c. 

The flux in the reflector has a different 
curvature from that in the core; but the 

core flux can be approximated by extending 

the point at which its characteristic function 
becomes zero by a distance such as ab, called 


the reflector savings.= 


REFLECTED 
REFLECT, 


"N-UNREFLECTED 






¢ 


r _— sass 

Figure 2 

Be Buckling, — p2, is a measure of the amount of 
curvature of the neutron flux pattem. The 
buckling which is required to make a given 
reactor critical can be found from the shape 
and dimensions of the reactor, and that which 
can be obtained by a given combination of 
materials can be found from the properties of 
the materials. The condition for criticality 


1. See Glasstone and Edlund, 8.25, and Appendix F. 
2. See Glasstone and Edlund, 7.20 to 7.60. 





of the reactor is that these two bucklings be 
equal. Buckling can also be used to find the 
proportion of neutron leakage from a given 
reactor. The probability of non-leakage dur- 
ing the slowing down process is e-BY and the 
ratio of thermal leakage to thermal absorption 
is L°B“, where L is the thermal diffusion 


ane 
length. 





Figure 3 


3. Thermal properties of the unit cell. 

a. Maximum heat generation rates qn: The rate 
at which heat can be generated in a fuel ele- 
ment is Limited by its heat transfer proper- 
ties. Along an element, the heat generation 
rate is proportional to the fission rate which 
is proportional to the thermal neutron flux, 
which will vary from zero at the ends of the 
reflector savings to a maximum halfway between 
them. For equal reflector savings at each 

- Glasstone and Edlund, 5.62 to 5.91 
See W. M. Rohsenow, J. Lewins , and J. P. Barger, "Steady 
State Temperature Distribution in a Nuclear Reactor with 


End and Genter Fed Coolant", NT=-269, MIT, hereafter cited 
as Rohsenow, Lewins, and Barger, pages 1-4. Also see Appendix G. 





end of the elements, this maximum will also 
be at the geometrical center of the active 
core. The coolant flows along the element 
and its temperature, t,, rises as it does 
so to a maximum at the outlet, while the 
temperature at the element centerline, t¢ , 
rises to a maximum at a point between the 
middle of the element's length and the out- 
let, then drops off due to the reduction in 
the heat generation rate. ‘he rate at which 
heat can be generated in a given element is 
limited by: 

(1) the maximum fuel centerline temperature 
must be maintained below the melting 
point of the fuel material, and, for 
uranium metal fuel, below the phase 
change at 660 °F. 

(2) the maximum coolant temperature must 
be less than the boiling point of the 
coolant. Boiling will change the water 
density, changing the density of hydro- 
gen atoms in the moderator and thereby 
the slowing down characteristics and 
the reactivity of the reactor. In ad- 
dition, boiling will reduce the heat 
transfer rate away from the element, 


leading to a temperature rise which may 


10. 
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be sufficient to melt it. 

(3) The temperature of the cladding will 
follow a pattern between those of the 
element centerline and the coolant. 

This must also be maintained below the 
melting point of the cladding. 
To improve the thermal efficiency of the plant, 
it is desirable to make the coolant outlet 
temperature as high as possible. Since this 
is limited by the boiling point of the coolant, 
the allowable temperature may be increased by 
increasing the water pressure. As can be seen 
from the curve of Figure b= preseures greater 
than 2000 psia give very little gain in tempera-— 
ture; moreover, the critical point of water, 
at which the distinction between the liquid 
and vapor phases is lost, is at 3206.2°F, which 
places an absolute upper limit on the pressure. 
A more practical limit is imposed by the dif- 
ficulties of manufacturing pressure vessels 
of the required sizes for pressures greater 


than around 2000 psia. 
Joo 





/000 | SLOOD no 
Prta 
Figure 4 
1. J. H. Keenan and F. G. Keyes, “Thermodynamic Properties of 
Steam", Wiley, (1936), page 39. 








lz. 


4, Thermal properties of the whole reactor. 


a. Overall heat generation rate.=° The heat genera- 
tion rate in each element will be proportional 
to its gq), which, in turn, is proportional to 
the flux distribution in a plane perpendicular 
to the axis of the element. Thus the maximum 
heat generation rate allowed for an element 
can exist only in one element, at the point 
of maximum flux, and all other elements will 
operate at heat generation rates less than the 
maximum. The ratio between the maximum and 
average heat generation rates and the maximum 
and average thermal neutron fluxes will be the 
same. 

b. Temperatures, t. The coolant inlet temperature 
is determined by the power system outside the 
reactor; the coolant outlet temperature is the 
temperature of the mixed outlets of all the 
channels. The average reactor and coolant 
temperatures can be found from the variations 
in point temperatures through the core. 

A uranium fuelled, light water cooled and moderated, pres- 
surized water reactor can be described by several basic para- 
meters: 

l. Describing the unit cell: 


a. Element shape and size. 


1. See Appendix G. 
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b. Volume ratio of moderator to fuel. 
ec. Material and thickness of cladding. 
d. Fuel material and degree of enrichment in 
U-235. 
2. Describing the reactor as a whole: 
a. Core shape and dimensions. 
b. Reflector dimensions. 
c. Temperature limitations. 
d. Power output. 

Changes in these parameters act andinteract in complex 
ways to affect the performance of the reactor. The number of 
parameters involved and the ranges of possible variation make 
experimental determination of the optimum combination for given 
requirements a long, tedious, and expensive process. 

Exponential experiment s=.. have been performed by H. Jd- 
Kouts at the Brookhaven National Laboratory on a number of 
the combinations*= and some of these have been continued to 
criticality by S. Krasik and A. Radkowsky (of the U. S. Atomic 
Energy Commission) at the Westinghouse Electric Company's Bettis 
Plant.2> Both these series involved uranium metal fuel elements; 
1. See Glasstone and Edlund, 9.90 to 9.111. 

2. Herbert Kouts, "Intracell Flux Traverses and Thermal Utiliza- 
tions for 1.15% Enriched Uranium Rods in Ordinary Water", 
Brookhaven National Laboratory Report No. 1987, August 11, 
1954. H. J. Kouts, J. Chernick, and I. Kaplan, "Exponential 
Experiments on Light Water Moderated 1 per cent U=-235 Lat= 
tices", Brookhaven National Laboratory Report No. 2094, 
November 28, 1952. 

3. S&S. Krasik and A. Radkowsky, “Pressurized Water Reactor 
Critical Experiments", paper no. 601, presented at the 


Geneva Conference on the Peaceful Uses of Atomic Energy, 


1955. 





@ similar, though shorter, series of critical experiments on 
uranium oxide elements has been made by J. R. Brown, A. Z. 
Krantz, and others at Bettist— 

These data are not by themselves sufficient to select 
an optimum configuration for a given set of requirements, but 
they are sufficient to provide experimental verification of 
calculations. 

The experimental results include critical bucklings, 
and, in those continued to criticality, critical number of 
rods, disadvantage factors, fast fission factors, resonance 
escape probabilities, and thermal utilizations for a series 
of rod diameters and volume ratios at a limited number of en- 
richments between 1 and 1.3%. In addition, Kouts! work in-= 
cludes sole extremely valuable plots of relative neutron flux 
along various lines through the cells. All the fuel was in 
rod form and clad with aluminum, with active core heights of 
48 inches. All cores were surrounded by effectively infinite 
water reflectors. 

For the purpose of theoretical investigation of a number 
of different configurations, several mathematical models are 
available, varying in both complexity and accuracy. The 
I. J. R. Brown, B. H. Noordoff, F. S. Frantz, J. J. Volpe, and 

D. R. Harris, "Measurement on Low Enrichments Light Water 
Moderated Critical Assemblies of U and UO2", paper no. 
XVII-5 presented at the Chicago Meeting of the American 
Nuclear Society, 1956. 

A. Z. Krantz, G. O. Smith, D. Klein, and W. Baer, "f, p, 
and Measurements on Slightly Enriched Light Water Moder- 


ated UO2 Lattices", paper no. XVII-6 presented at the 
Chicago Meeting of the American Nuclear Society, 1956. 
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simplest of these, the “one group" squatiaee assumes that all 
neutrons are at thermal energy. The critical condition is 
Keo 


1 ft B“Le 


This is relatively simple mathematically, but poor in accuracy. 


Somewhat better in accuracy and not much more difficult in 
computation is the “modified one group" theory which in= 
cludes a factor for leakage during the slowing down process. 


The critical condition may be expressed as 


— » ~BeT 


= 2 22 
lL + BL 


the Fermi Age variant, or 


E: Ke 
i + 

1 t Bem 
where M*, the migration area 2° is (Le +7), for a large 
reactor. 

a2 
More complex is the “two group" theory; which assumes 

a fast and a thermal group of neutrons with average reactor 


properties for each. Four differential equations are formed, 


expressing fast and thermal flux in the core and the reflector 


1. See Glasstone and Edlund, 8.15 to 8.37. 

2. See Glasstone and Edlund, 7.20 to 7.27, and 7.62 to 7.64 
3. See Glasstone and Edlund, 7.63 , 

4. See Glasstone and Edlund, 8.38 to 8.58 


is5.c 


De W Or : Z10c 1 i kate, P2, 
Dog VW Pac ~ = 22/0 be t Fre Vie = 
Dip 7 Orr i Zir Dir - 
Dar VW Per - Zar Por ‘ Zi Qir = 


Where the subscripts 1 and 2 refer to the fast and thermal 
fluxes respectively and the subscripts c and r refer to the 
core and reflector respectively. 

These are solved by means of boundary conditions which 
require both fluxes to be continuous across the core-reflector 
boundaries and both fluxes to go to zero at the outer reflec- 


tor extrapolated boundary. These solutions lead to a critical 


determinant, 
x Y -2) 0 
51x g5¥ -$32, “Fe 
Pie *" Dy oY! Pir4y of 
5 D2,%" B2D2¢%! -83D2r2y | “D222 


where, for a circular cylinder with an infinite radial reflec= 


tor, 
x = Jo( Mr). y = Io Wr) 
5) = _ x = ae, So = Dic © Pp 
Do 2 D 1 Z 
Cue’ era Nov2e y= V 
L2e L2¢ 
2 ls MylKi vr) = Grr 
2 oN 3 D 2 2 
ma (% oy - Wy 
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p = diffusion coefficientt 

Tf ts Fermi Age 

Coefficients for other geometries are given in Glasstone and 
Edlund, sections 8.46 and 8.51. 

This treatment is capable of considering only one reflec- 
tor. But in practical reactors, reflectors are usually pro- 
‘vided for all exposed faces. This situation can be handled 
with moderate accuracy by Hill's Approximation=~ in which the 
reactor is considered separately for each reflected dimension 
and the solutions superimposed. In this approximation, the 
corners are not properly considered, but the flux is generally 
low in those areas and the resulting error not excessive. 

It is possible to reduce the labor involved in critical 
computations considerably by use of a digital computer. These 

calculations are particularly useful in an investigation such 
as this, where a large number of similar problems will have 
to be solved by identical methods, for once the program is 
written, it cari be used to solve any number of cases. 
“‘Bwo such programs for Whirlwind, the high speed digital 
computer at MIT were available, both written by James R. 


Powell as part of his MIT thesis for the Sc.D. degree. = 


1. Glasstone and Edlund 5.7 and 5.20 to 5.24 

2. Powell, James R., “Nuclear Characteristics of Heavy Water 
Moderated, Enriched, Homogeneous Reactors”, Sc.D. thesis, 
Mans L957 





These programs were based on the assumptions of a heavy water 
moderator, highly enriched uranium, and a homogeneous core. 
One of these programs was a computer representation of the two 
group Hill‘’s approximation method already discussed; the other 
was based on a Fourier treatment of the diffusion equations. 
The theory of the Fourier treatment will be briefly out- 
linedhere.— A detailed description of the theory for the 
Fourier program is contained in Powell's thesis. The thermal 
diffusion equation states that the losses due to leakage and 
absorption must equal the number of neutrons s lowing down to 


thermal energy for a critical reactor. 


va - Seine! 8 = 0 (1) 


Leakage absorption production 


Each of the terms in equation (1) can be expanded in a Four-= 
ler series with eigenfunctions obtained from the wave equa- 


tion 

7x, 4 Be, 0 
to satisfy boundary — on: This thermal flux expansion 
would take the form of ; 


1. Based on J. R. Powell's Thesis, Section III. 
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Since absorption and slowing down vary across the reactor and 


depend on factors other than flux, two further expansions are 


ko (r) Zalr r - 
Keo (2) Zale) Pir) Z eyy(r) 


Since &,(r) is not constant, but varies as a step function 


required 


of (r) from region to region an additional equation is write 


ten relating the &, with the regions. The coefficients of 
these expansions are then related to the Fourier coeffici- 
ents of the thermal flux expansion. Using these relations, 
an infinite series of equations is obtained with one equation 
for each unknown Fourier coefficient. For non-trivial solu= 
tions the determinant of the coefficients must equal zero. 
This constitutes the critical equation. In the program the 
critical mass is found using a limiting process, by finding 
the mass of U-235 needed to make the determinant of the first 


two rows and columns vanish, then the determinant of the first 


19. 


three rows and columns vanish and noting that the mass approaches 


a@ limiting value as the number of rows and columns increases. 
Powell points out two major advantages to the Fourier 
solution of diffusion equations. The first of these is the 
ability to represent the slowing down behavior of fast neu- 
trons as a convolution of an age and diffusion kernel. The 


Fermi age kernel is expressed by 


Pj = exp( =A net) 
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and the diffusion kernel is expressed by 


Po = a eT 
Ie J, 8,° 
By postulating that the neutrons followed Fermi age theory 
until they reached energy Ey and then diffused away from the 
source until a collision reduced their energy to thermal 
energy, E,y,, it was felt by J. Powell that effectively two 
fast groups of neutrons were being considered. The slowing 


down kernel was then represented as 
Z | 2 
Ry [ tm, By). Ba | Po (Ey By ) 
a Qepeee 
‘ i liga 
1 +t “ 
We By 


The relative size of Fy and 2 was determined by matching 


the observed distribution of slowing down density with the 
calculated distribution so that no definite energy boundary 
between fast groups is established. 

The second major advantage of the Fourier treatment 
pointed out by Powell was the relative ease with which addi- 
tional reflectors could be handled. As a result of this, 
the Powell program included the provision for two top and 
two bottom reflectors and for six radial regions, any one 
of which could be the core region. All cores in this thesis 
are analyzed or designed on the basis of two regions only, a 
core and an outer infinite water reflector on top, bottom and 


sides. 





A choice of two methods of machine computation was thus 
offered: the Fourier treatment and the two group Hill's ap-= 
proximation. For several reasons, the Fourier program seemed 
preferable for the purposes of the thesis, and it was ulti- 
mately selected. These reasons were: 

1. The simultaneous solution of two cases required 
for Hill's approximation makes the solution 
rather lengthy. Computer time required for this 
program is about four minutes compared to two 
minutes for the Fourier. 

2. Comparisons made by Powell in his thesis of the 
calculated results by both methods against ex- 
perimental data for highly enriched, heavy 
water moderated cores showed smaller errors in 
the critical masses calculated by the Fourier 


program. Errors in cylindrical geometry varied 
from t -08 to - 14% with the Fourier and from 


~ 4.0 to - 20% with the two group Hill's approxi- 
mation. 

3. Use of a two group program requires inputs of 
fast group constants. These are not well known, 
and, since a variety of materials was to be 
used in the core, it appeared that the introduce 
tion of this added uncertainty would not be ad= 
vantageous. The Fourier program, on the other 
hand, uses only thermal constants and a Fermi 


Age, all of which are known to fair accuracy. 


Ale 
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Having the required program inputs, the data available 


for experimental checks, the practical problems of fabrica— 


tion which would limit the choices in a real systen in mind, 


the following ranges of parameter variation were selected: 


Bh 


In the unit cell. 

a. Fuel elements to be rods, from 0.1 to 0.6 
inches in diameter. 

bd. Volume ratios, moderator to fuel, from 1.5 to 
6.0. 

c. Cladding to be .015 inches thick, both stain- 
less steel and zirconium to be considered. 

dad. Fuel to be uranium oxide, UOo; quantity and 
enrichment to be determined. 

For the whole reactor. 

a. Cores to be right circular cylinders not to 
exceed 13 feet in diameter, and to vary between 
4 and 8 feet in height. 

b. Cores to be surrounded by an 18 centimeter thick 
water reflector both radial and top and bottom; 
this is effectively an infinite reflector. 

ce. Maximum temperature in fuel: 4800°F 
Coolant pressure: 2000 psia 
Coolant temperatures: Maximum: 628°F 

Inlet : 487°F 
Average: 508°F 


d. Thermal power output: 400 megawatts. 
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In order to carry out the parametric study utilizing 


the machine calculation, a series of steps were undertaken: 


1. 


A method was developed to convert the program from 
the high enrichment, homogeneous case with heavy 
water moderator and aluminum structure for which 
it was written, to the case which was to be in- 
vestigated: low enrichment, heterogeneous, 

light water moderated, and stainless steel or 
zirconium structure. 

A method was developed for designing a consistent 
series of cores for testing which reflects the 
influence of parameters to be varied both on the 
heat transfer and the nuclear characteristics 

of the core, and supplies the necessary data for 
entering the machine calculation. 

A survey was made to collect best values of the 
necessary constants and basic data, and to select 
methods for the various preparatory calculations 
required. 

Correlations were obtained, both for the prepara= 
tory calculations and for the computer results, 


with the experimental data available. 





il. 
PROCEDURE 


The use of the Whirlwind program selected, which was 
written for a highly enriched, heavy water moderated, homogen- 
eous reactor, for computations on a low enrichment, light water 
moderated, heterogeneous reactor involved a considerable amount 
of calculation prior to entering the machine computation. 

The program assumes the existence of only three materials 
in the core: 

1. Fuel, assumed to be U-235. 

2. Moderator, assumed to be Do0. 

3. An absorber, assumed to be aluminum. 
Appropriate values for the densities and atomic weights of 
these materials are built into the program. Further, the 
program assumes that the core region consists of a uniform 
mixture of moderator and absorber with a small amount of fuel 
uniformly distributed in the mixture. Among the computer 
inputs? - are values of the microscopic absorption and transport 
cross sections of moderator and absorber, the volume ratio of 
absorber to moderator, and the microscopic absorption and 
fission cross sections of the fuel. The principal problem 
in utilization of the program was to compensate within these 
inputs for the effects of heterogeneity, the use of materials 
other than heavy water and aluminum as moderator and absorber, 
and the additional presence of relatively large quantities of 


U-238, for which the program makes no provision. 


I. See Appendix &E. 
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In principle, two reactors, one heterogeneous, the other 
homogeneous are equivalent if the average values of macroscopic 
absorption, fission, and transport cross sections for the en- 
tire reactors are equal. This is the basis of the derivations 
for the necessary adjustments to the cross sections as outlined 
in Appendix B. Each value of a thermal cross section is ad- 
justed to have the same value as the effective volume average value 
of the corresponding cross section in the heterogeneous reactor 
weighted by the relative flux through the unit cell. For ab- 


sorption in the fuel, as an example, this is expressed as 


follows:? 
S (2 5) = N(25)0-° (25) Yo® «x N'(25)o~a! (25) (1) 
Yo gd, +t Vi 
E a(2g) = Nl25)eg(25) VO = NI25) Tere PW (2) 
o-ma'(25)= oa(25) Von : (3) 
VotVi G 


o 
where N(25) and o~(25) are the number of atoms per cubic centi- 


meter of fuel and the microscopic cross section in the hetero- 
geneous reactor and N'(25) ando-4'!(25) are the number of atoms 
per cubic centimeter of core and the microscopic cross section 
in the equivalent homogeneous reactor, g-4'(25') is the appro- 
priate computer input. 

Af 4s the mien factor, discussed in the 
Introduction. Methods of computing the disadvantage factor 
are discussed in Appendix A. The choice of method for use in 
this thesis was based on the results of A. L. Kaplan's M.I.T. 


I. See Appendix B. 
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thesis, *- in which he compared the accuracy of the disadvan- 
tage factor for the M.I.T. Reactor as computed from three in- 
creasingly sophisticated models. His conclusion was that the 
most complex of the three models offered only a three per cent 
increase in accuracy over the simplest, and that the added 
accuracy did not justify the added labor of computation. The 
model chosen, then, was the simplest of Kaplan’s three; it 
assumes constant fast flux across the cell and no net thermal 
leakage from the cell.*: 
Since the program is based on pure U-235 as the fuel, 
it appeared that the best way to include the effects of U-238 
was to consider it as part of the absorber, which then consisted 
of both U-238 and the cladding material. The process developed 
for adjusting the constants for this arrangement was as follows: 
1. A macroscopic cross section was formed for the 
U-238 from N(28) in the fuel andoq!(28), obtained 
by a disadvantage factor adjustment as in equation 
(3) above. 
2. A macroscopic cross section was obtained for the 
cladding material. Since the cladding is in a 
region of the unit cell in which the flux is very 
nearly the average value, no disadvantage factor 
correction was included, and the macroscopic 
cross section consisted simply of the product 


N(clad) xc~4(clad). 


1. A. L. Kaplan, “Theoretical S+udies of Neutron Flux Distri- 
bution Expected in the MIT Nuclear Research Reactor", MIT 
S.M. Thesis, 1955. 

2. Also discussed in Glasstone and Edlund, 9.47 to 9.66, and 
in Appendix A. 
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3. These cross sections were weighted by their respec= 
tive volumes and averaged over the total volume 
of fuel and cladding. 

4. Since the machine would compute a macroscopic ab 
sorber cross section on the basis of the density 
and atomic weight of aluminum, which were integral 
parts of the program, the average macroscopic 
cross section obtained in step (3) was divided 
by N(Al). 

Mathematically, the process is described by the following equa- 


tion: 
* 2 '(2 d 
c(t absorber) N(28) (28) xe Sa eee V(clad) , 


Aal 


fe ly 
fa x 6.025 x 10£3 am 


A step for spreading the absorber uniformly through the core 
volume is part of the computer program; a program input, "x", 
is the volume ratio of absorber to moderator. For computing 
"x", the volume of fuel rod plus cladding was used. 

A further adjustment must be made for the presence of 
U-238 in the fuel rods. With highly enriched uranium, there 
is neither a fast fission effect nor resonance absorption, and 
the infinite multiplication factor, k,, is simply 

Koo = VHF 

Yfor the fuel 1s a program input; f 1s computed by the machine. 


However, for low enrichments, 


| (apa 











28. 


In order to use the program for low enrichments, then, € and 
p were separately calculated, as described in Appendix D, and 
the product was used as the program input "7". 

An additional arbitrary factor of 1.1 was used in the 
program input 17m to bring computed results into agreement 
with experimental data. The derivation of this factor is dis- 
cussed in Appendix H. Thus, 

"9a" =7(25)ep x1.l 

Similar adjustments were made for the substitution of 
light water at 508°F for the heavy water at room temperature 
on which the program is based. Here the adjustment process 


was as follows: 


0 V1 + Vo A ps0 
g7(moderator) = H20 5» O- #20 4 ae » els 
AH20 Vit Vo “4 
| 1 
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Appendix F lists all cross sections required to calculate 
both the cores for comparison with experiment and cores for 
our parametric study. Cross section information was obtained 
from BNL=-325, Neutron Cross Sections, and modified for Maxwell- 
Boltzmann distribution and additionally for temperature where 
necessary. As pointed out in the Introduction, one of the ad- 
vantages of the Fourier program is its ability to use both a 
Fermi age and a diffusion age and in this manner effectively 
treat the fast neutrons as two groups. A rather complete sur- 
vey of the literature was made to determine whether any theoreti- 
cal work had been done on defining age in this way for light 
water, and it appeared that it had not. The most authoritative 


experimental results were computed by Kouts, et.al. at Brookhaven 





based on the Fermi Age theory. These results indicated that 
7 = 31 em@ for light water of density 1 gm/cm3 for water to 
metal ratios from one to three. It was therefore decided to 
use this number as the Fermi age and to set the diffusion age 
input equal to zero thus eliminating this. The resulting cal- 
culation performed by the program is actually a one group con- 
tinuous slowing down solution. Appendix F inoludes a fuller 
description of this literature survey. 

In addition to the neutron characteristics of the mate-~ 
rial components of the reactor, Powell's program also requires 
the dimensions of the reactor as inputs. Since uniform reflec- 
tors and cylindrical geometry were assumed for all cases, the 
only dimensions required were core height and radius. A series 
of set core heights was assumed; for given height, power out- 
put, and cell geometry, the radius is determined by the heat 
transfer characteristics and temperature limitations of the 
materials and the overall flux distribution. 

A method of computing core radius on the basis of cell 
geometry and materials for given power output and core height 
was developed from the functions derived by Rohsenow, Lewins, 
and Barger? ° for temperature distributions in a core. This 
development, which is carried through in detail in Appendix G, 
assumes that the optimum configuration is that in which the 
temperature limitations for both bulk boiling in the water and 
melting of the fuel element are just met at the hottest points 
in these media. A third possible criterion might have been 
that the melting point of the cladding should just be reached; 


1. Rohsenow, Lewins, and Barger, op.cit., and Appendix G 
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the temperature of the hottest point in the cladding was 
checked in each calculation and found to be well below the 
melting temperature; thus this was not a limiting factor. 

This method offers considerable advantages ina study 
such as this. By using it, it 1s possible to design a series 
of cores of widely varying geometry all of which will be 
capable of developing the same amount of thermal power while 
just meeting an identical set of temperature limitations. 
This allows meaningful comparisons between the different com- 
dinations of parameters. 

The process carried through on each of the cores con- 
sidered was as follows: 

1. Characteristics were selected from the ranges of 
parameter variations listed in the Introduction. 
These were: rod diameter, moderator to fuel 
volume ratio, cladding material (stainless steel 
or zirconium), and core height. 

2. Using these characteristics, adjusted cross sec- 
tions were computed. 

3. Using the same characteristics, a core radius was 
computed; since this was the boundary of a uni- 
form lattice of cells each of known geometry, 
the number of cells and hence the total amount 
of fuel was also known. 

4, Using the results of steps (2) and (3), as input 
data, a Whirlwind computation using Powell's 


Fourier program was made to find the critical 
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mass of U-235 required. 

5. The critical mass of U=-235 was compared to the 
total mass of uranium, found in step (3), giving 
the required enrichment of the uranium for critical- 
ity in that core. 

6. Cores were compared on the basis of amount and 
enrichment of uranium required. 

It will be noted that enrichment is required in order 
to find the disadvantage factor, and therefore in order to 
compute the adjusted cross sections; it is also required to 
determine the amount of U=-238 present and hence the U-238 con-~ 
tribution to the "absorber" cross sections. A solution can 
thus only be obtained by an iterative process in which an en- 
richment 18 assumed in step (2) and checked in step (5). 
However, since the range of enrichments involved is very 
small, from about one to three percent, and the solution 
converges very rapidly, only a few cores had to be run more 


than once, and none more than twice. 
























TABLE I 
Total Mass of Uranium (106 om ) 
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TABLE IT 


Stainless Steel Clad 
U-235 Critical Mass (10+ gm) 
Diameter] Height 
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TABLE III 


Zirconium Clad 
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TABLE IV 
Enrichment 
Stainless Steel Clad 
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TABLE V 
Enrichment 
Zirconium Clad 
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TABLE VI 





Core Diameter to Height Ratio 
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TABLE VII 
ko 


Stainless Steel 
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pti ol RESULTS 
The results of the calculations for the series of 400 

megawatt cores described in the Introduction are shown in 
the Results, and total uranium mass and enrichment are plotted 
in Figures 2, 3 and 4 Both the total mass of uranium and 
the required critical enrichment are shown, since the uranium 
cost is a function of both. | 

Total mass of uranium is a function solely of the heat 
transfer properties of the combination of parameters under 
consideration. It is a strong function of rod diameter and 
a weak function of core height and moderator to fuel volume 
ratio. For a given core height the allowable heat generation 
rate in the central rod is almost independent of rod diameter. 
This is due to the low thermal conductivity of the oxide 
fuel (0.95 Btu/hr.-ft.-°F) compared to those of stainless 
steel (10.6) or zirconium (8.15). As a result, almost the 
entire temperature drop between the 4800°F fuel centerline 
and the water at less than 628°F is taken across the fuel 
region. Using the methods of Rohsenow, Lewins, and Barger, 
a sample temperature distribution was computed for the 
point of maximum fuel centerline temperature in a 0.4 inch 
diameter rod. It is shown in Figure 1, in which the magni- 
tude of the drop in the fuel, in this case 97% of the total 
drop, can be seen and compared with those in the clad and 


across the filn. 


1. Rohsenow, Lewins, and Barger, op. cit. 
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Figure 1 


| 1. 
Assuming that the flux in the fuel can be expressed as: 


? = alo(Mr) (1) 
then the heat generation pattern will be the same: 


a(r) = AI,(K or) (2). 
where q(r) is the volume rate of heat generation at radius r. 
The average heat generation rate for the whole rod, q, can be 


expressed as: 


‘fn 27rr Aly ( #or)ar 
0 


0 


TH Ro 
211 (Mo Ro) I, (6 r) (3) 


z 
TTRo 


and q(r) 


The total heat generation Q(r), inside any radius, r, will be: 


1. Glasstone and Edlund, 9.57 





Figure 2 
Total Mass of Uranium 
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Figure 4 
Enrichment vs Moderator to Fuel Ratio 
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> i 
a) = Sen fe Be I a 
| , 27; (2 Ro) o( Ko rjdr 
= 714 FR, ri, ( r) (4) 
TiVKoRy) 1 (Fo 


For steady state conditions, all this heat must be transferred 


across the circumference at r. Thus 


Q(r) = Th Ro = wt. ( Aan) = - k.2nr.dt/ar 
; Th (XK oRo) 1 Ao 





at/ar = Fo i - Ty (Xor) (5) 


and integrating: 


to -t(R,) = At = BRD) is * [Fol HoRo)Y (6) 


Since the total heat generated, Q, is To Ro” , 


1 


_ Ro se (me Roel 
at = Ro is - 11 ( Te Ro) 
ar = 2 x a * So)? (7) 





aT Ti ( Xo Ro) 


The third factor in equation (7) is very nearly constant over 
the range of @>Ro involved, ranging from 0.499 at HRo = 0.05 
to 0.459 at QWRo = 0.40. Thus, since 271Tk is constant, for 
constant At, Q is the same regardless of the value of Rp, 
and very nearly the same power per unit length is developed 
in the central rod regardless of its diameter. | 
Figure 5 shows the effect of increasing the radius of 


the core for the same central rod heat generation, assuming 
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CORE RADIVS 
Figure 5 


a J, overall flux distribution. The average to maximum heat 
generation ratio is reduced slightly, due to the large nunm- 
ber of rods near the outer circumference which operate at 
low powers. This difference 1s, however, small; the ratios 
vary between about 0.5 and 0.6. 
The two influences that will tend to increase the core 
radius for constant power output are: | 
1. {Increasing the volume ratio of moderator to 
fuel. 
2 Decreasing the core height. 
The result is the pattern observed in Figure 2 and in 
Table I of the Results: since very nearly the same power per 
unit length can be developed in the average rod regardless 
of rod diameter, very nearly the same number of rods is re- 


quired at a given core height regardless of rod diameter, 
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and the mass of uranium increases roughly as the square of 
the rod diameter. In the range of core parameters investi- 
gated, which represents about the range of practicality, the 
number of rods for given core height varied through about 
20%. And since increase of the moderator-fuel ratio inoreases 
the diameter, a slightly Larger number of rods is required, 
and, therefore, the mass of uranium is slightly increased. 
Similarly, a reduction in core height demands an inorease 
in the number of the shorter rods, resulting in a diameter 
increase, and, again, a lowering of the average to maximum 
heat generation ratio and a slight increase in the mass of 
uranium. 
The variation of enrichment with rod diameter and 
moderator to fuel ratio is a somewhat more complex phenomenon. 
The effect of leakage is very small, for two reasons: 

1. Both Fermi age,F? , and the squared diffusion 
length, Lé, are small for water moderated lat- 
tices. Room temperature values are 31 om and — 
4 om respectively; both increase with a decrease 
in water density, so that at the test conditions, 
the values are 49.5 em@ and 5 cm*. 

2. Due to the heat transfer requirements, the cores, 
of whatever proportion of length to diameter, 
are relatively large and the buckling consequent- 
ly small. For the cores considered, B* computed 


from the core dimensions bye 


1. Glasstone and Edlund, 7.60 
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2 2 
2. {2-405 + 
_ ( R ) 4 
varied between 3 x 10-4 ana 9x 10-4 
From these values the non-leakage probability can be com- 


puted approximately from the modified one group theory dis- 
cussed in the Introduction, as the product of slowing down 





: and thermal diffusion non-leakages: 





~BeY 


6 


slowing down non-leakage 


1/(1 + L@B@) 


total non~leakage = e~B2P / (1 + L°B~) (8) 


thermal diffusion non-Leakage 


For the range of bucklings and the appropriate age and 
diffusion length values, the non-leakage probability computed 
in this way for the cores considered varies between 0.95 
and 0.985. This is so small a range as to be nearly constant. 

The machine calculation gives as a supplementary result the 
value of Koo required for criticality; these are shown in 
Tables VII and VIII of the Results; it will be seen that they 
vary, for stainless steel cladding, from about 1.06 to about 
1.13. Since 

k (effective) = 1 = Kyg x non-leakage (9) 

the non-leakage computed by this somewhat more sophisticated 
method varies in the narrow range from 0.88 to 0.94 still very 
nearly constant. 

This means that leakage alone will not explain the 
strong variations found in enrichment, which must then depend 
on the effects of the cell parameters. these effects can be 


explained qualitatively with the use of a somewhat simplified 
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model, as described below. 
As discussed in the Introduction, kyo is the product 

of four factors; a fifth was introduced to bring the computer 
results into agreement with experiment ,*— so that 

koo = 1.17€pf 
Of these factors 7 is a constant for U-235 fuel 

kee = l-lex 2.08 ¢ pf 
Of the remaining three factors, € and p were computed in the 

2 


work preliminary to entering the machine calculation=" and 


f was, in effect, computed by the machine. fhus, 


k - 
ce mel sg Fee 2 = 56 Tee (10) 
Llex 2.08 Ep Ep 


By equation (8), 
Koo = 1/(total non~-leakage) 


And, using an average value of 0.90 for non-leakage, 


f= 0.436 x 5G x it = a pee (11) 


The thermal utilization, f, is also definied as the 
ratio of absorption in the fuel to total absorption in fuel, 


U-238, and moderator: 


Yo Sa25 Po (12) 


Vo a25 Yo + 4, = a28 Vo t Fy 201 Pi 


Algebraic manipulation of this gives 


1 = Yo Sa25Go 4 Vo Za28do + Vv, Eal f1 
: Yo Za25 Po Vo &a25 Qo Vo Fa2s5 Po 


{ 2 


Li i} 28,028 + 2 pie op Sige (13) 
r N25 025 Vo Po FS a25 


lL. See Appendix H. 
2. See Procedure and Appendix D. 








For low enrichments, 
enrichment = Nes/(Nos5 + Nog) ~ Nos/Nog 
And for enrichments in a narrow range around 2% and cladding 


volume much smaller than water volume, 


TEE /ag, = constant ~<~ 0.1 
Substituting into equation (13) 

1-2 + oc x 2 xGVil + 4.15 x 1073 

r - enrichment 
And from equation (11) 


Lé 
¢> oye, = 2-06 ep 


Equation (15) may be plotted, as in Figure 6a from the com- 


bl. 


(14) 


(15) 


puted data on fast fission factor and resonance escape proba- 


bility, which are functions only of rod size and moderator 
to fuel ratio. } 

The two terms of equation (14) can be separated; the 
first term is a straight line of slope 0.1 with increasing 
volume ratio, slightly modified by the disadvantage factor, 


QP, 6 » which was also computed in the preliminary calcu- 


lations for each core and increases slowly with increasing 
rod diameter and volume ratio. this is also plotted in 
Figure 6a. The difference between the two lines so plotted 
for each rod size must be the second term of equation (14), 
and therefore proportional to the reciprocal of enrichment. 
It will be seen from Figure 6b, in which these differences 
are plotted for the rod diameters against the moderator to 





Figure 6 


Effect of Cell Parameters on 
Enrichment (Equations (14) and (15)) 
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Moderator to Fuel Ratio , V1/Vo 
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fuel ratio, that the difference reaches a maximum, represent- 
ing a minimum enrichment, for each rod diameter, and that 
these maxima are reached at a different moderator to fuel 
ratio for each rod diameter, which increases with decreasing 
rod diameter. 

The same pattern can be seen in the results of the more 
accurate machine computation, as shown in Figures 3 and 4 and 
in Tables IV and V of the Results. Optimum enrichments are 
found for each rod diameter at moderator to fuel ratios which 
increase with a decrease in rod diameter. If leakage is as- 
sumed small and constant, the enrichment depends on the cell 
effects shown in Figure 6a. The product ¢p is greater for 
larger rod diameters, and increases with decreasing slope 
as the moderator to fuel ratio increases. However, since 
water is a strong absorber of neutrons, the absorption effect 
increases steadily as more water is added, tending to resist 
the advantage of increased ¢p. In addition, as rod size in- 
creases, the disadvantage factor effect becomes larger, also 
tending to resist the effect of €p. As a result, an opti- 
mum moderator to fuel ratio is reached, beyond which a fur- 
ther increase will result in an increased enrichment. Due 
to the earlier leveling off of ¢p and the greater effect of 
disadvantage factor, both of which can be clearly seen in 
Figure 6a, this optimum occurs at lower values of the volume 
ratio for larger rods. 

That this effect is nearly independent of core geometry 
at this power is shown in Figures 7 and 8a in which enrichments 
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CORE DIAMETER TO HEIGHT RATIO 
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are plotted for rod diameters and core heights against the 
diameter to height ratio of the core. If core geometry had 
a significant effect, it would be shown by the existence of 
a minimum enrichment near a diameter to length ratio of one. 
Although slight variations exist in the best enrichment ob- 
tainable with a given rod diameter, they are significant 
mainly in that they are so small. As indicated by Figures 
7, 8a and 8b, an optimum enrichment can be obtained for a 
given rod diameter which is independent of the proportions 
of the core over a wide range. 

Although the volume of cladding is small compared to 
those of fuel and water, its significance is seen in the conm- 
parison between the stainless steel and the zirconium clad 
fuels, shown in Figures 3, 4,7 and 8. Due to the much smal- 


ler neutron absorption in zirconium ( 2, = 0.00526 compared 


to | a = 0.1535 for stainless steel), criticality can be 


attained with considerably lower enrichments. That the. 
optimum volume ratios are not far removed from those found | 
for the stainless steel clad cores is due to the lowering of 
 Sa25 with the decrease in enrichment from about 2% to about 
1%. Thus = .3/ F225, the slope of the lower group of lines 


in Figure 6b before the disadvantage factor correction is ap=- 
plied, is not greatly changed, and the optimum enrichments 
occur at roughly the same values of Vi/Vp.- 

It should be observed that the minimum enrichment does 


not necessarily indicate the optimum core. The uranium cost 
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is a function both of enrichment and the total quantity of 
uranium. Since the total quantity of uranium goes up rapidly 
with the rod diameter and the enrichment is reduced only 
slowly, Atomic Energy Commission price figures released in 
December of 19562: indicate that the increased amount of 
uranium required for the larger rods will overweight the low- 
er unit cost due to the lower enrichment, making the fuel 
cost less for the 0.2 inch rods than for the 0.6 inch. This 
effect will be reinforced by the greater amount of cladding 
material required to place an equal thickness of clad around | 
the larger rods, since the number of rods is roughly the same. 
The cost of rod fabrication and of the pressure vessel are 
also significant items. 

Finally, the enrichments computed are those required 
to reach initial criticality only. In order to compensate 
for the burnup effects of operation at power over a long period 
of time, these enrichments would have to be increased by a 
factor which depends on the lattice parameters. A method of 
finding this factor is derived in Appendix J. The cores com- 
puted in the thesis would be the basis of such a calculation. 


LT ts! 


1. "Enriched Uranium Prices", Nucleonics 14, 12, R-2 (Dec. 1956) 


och eneiieente 
I. Relating to methods of calculation. 

1. Minimum core size can be computed for a specified 
power output and specified unit cell parameters. 

2. The Whirlwind program written by J. R. Powell using 
a Fourier method of solution for critical mass of a 
cylindrical, highly enriched, homogeneous, heavy 
water moderated reactor can be adapted for use in the 
cylindrical, low enrichment, heterogeneous, light 
water case, and experimental correlations were ob- 
tained for the adjusted program. 

3. The factor by which critical enrichment of a particu- 
lar core must be increased to attain a specified power 
time of operation can be computed as a function of the 
critical amount of U=-235 and the cell parameters. 

II. Relating to the series of 400 megawatt thermal output 
cylindrical pressurized light water reactor cores studied. 
1. Total mass of uranium required is a function of the 

heat transfer characteristics of the unit cell; it 
is approximately proportional to the square of the 
rod diameter. 

2. Enrichment required for initial criticality is princi- 
pally a function of the unit cell parameters and is 
nearly independent of the core size and shape, between 
diameter to length of core ratios of 0.8 and 3. 


3. For each rod diameter and cladding material an optimum 
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moderator to fuel volume ratio exists at which the 
eritical enrichment is a minimum. For the diameters 


tested, these were as follows; 








Rod Diameter Optimum V1/Vo 
Stainless Steel clad Zirconium clad 
0.6 inch 3.0 3.0 
0.4 inch 3.75 4.0 
0.2 inch We 5 e5 
An optimum rod diameter at which minimum critical en- 
richment is lowest exists for rods clad with each 
cladding material. For zirconium clad rods, this oc- 
curs at about 0.45 inches rod diameter*for stainless 
steel clad, at slightly over 0.6 inches. 
Rod Diameter Minimum critical enrichment 
Stainless Steel clad Zirconium clad 
0.6 inch 99 890 
0.4 inch Pole 085 
0.2 inch 1.5 Leis 
Critical enrichments are lower for zirconium clad rods 
than for stainless steel clad, of the order of 1% en- 
richment compared to the order of 2%. 
Optimum critical enrichment does not necessarily indi-~ 


cate the optimum core; other factors which must be 
considered are total quantity of uranium, cost data 
for uranium, cladding materials, and fabrication of 
fuel elements and pressure vessels, and the additional 


enrichment necessary to attain a specified core life. 
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APPENDIX A 
DISADVANTAGE FACTOR 


A. L. Kaplan in a Massachusetts Institute of Technology 
$.M. thesis titled "Theoretical Studies of Neutron Flux Dis- 
tribution Expected in the MIT Nuclear Research Reactor" com- 
pares the disadvantage factor calculated by two group diffu- 
sion theory based upon three reactor models. The first model 
assumes fast neutron flux constant through the cell and no 
net leakage of thermal neutrons across the outer boundary of 
the cell. This is identical to the development of Glasstone 
and Edlund. Model number two allows the fast neutron flux 
to vary throughout the cell as a function of the radial dis- 
tance from the center of the fuel element and again assumes 
no net leakage of thermal neutrons across the outer cell 
boundary. The last model differs from the second by assuming 
a net leakage of thermal neutrons such that the ratio of 
thermal leakage to fast leakage equals the ratio for an equi- 
valent homogenous cell. He concluded that the third model 
gave only three percent increase in accuracy over the first 
model and that this added accuracy did not justify the in- 
creased complexity in the calculation. This thesis work by 
Kaplan on heavy water, highly enriched cores and the close 
agreement between calculated and experimental disadvantage 
factors for light water, low enrichment cores (Table A-I) 
indicated that the one group diffusion theory in Glasstone 


and Edlund was adequate for this thesis. 
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This development approximates the reactor core by divid- 
ing it into a number of identical unit cells so that a square 
cross section is replaced by a circular cross section of equal 


area with radius R,; in the manner shown in the following figure. 





The equivalent cell then consists of a central fuel rod of 
radius Ry with a concentric volume of moderator (water and 


clad) as shown below. 











The Glasstone and Edlund development assumes slowing 
down density constant in the moderator and zero in the uraniun. 
With the boundary conditions of continuity of neutron flux and 
neutron current density at the fuel rod-moderator boundary 
and no leakage at the outer cell boundary, the following equa- 


tions for thermal utilization are developed: 


1+ yest Pe (e- 1) (1) 
r= Xo Ro Io (oRo) (2) 
Ii (HoRo) 





Eliminating thermal utilization from equations (1) and (4) 





: = = cD Vo 


where ¢, and “) are the volume average thermal fluxes in the 
fuel and moderator, respectively, and their ratio is the dis- 
advantage factor. This appendix contains a sample calcula- 
tion for this factor. 

In order to determine the accuracy of this method of 
calculation, a series of experimental flux traverses for 
light water moderated, low enrichment uranium lattices ob- 
tained at Brookhaven National Laboratory by H. J. Kouts, et. 
al. was utilized for comparison. Tables A-II, A-IITI and A-IV 
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and Figure A-I are excerpts from these reports. Table A-I 
compares calculated and experimental disadvantage factors 

over the range of experimental flux traverses presently avail- 
able. It should be noted that this is an extremely narrow 
range; three enrichments, 1.027%, 1.15% and 1.30% for 0.6 
inch rods and one enrichment, 1.027%, for 0.75 inch rods. 

Agreement between calculated and measured disadvantage 
factors for the 0.6 tnch rods are all less than 6% error. 

With increased diameter the percent error becomes greater. 
This is explained by the failure of the theoretical functions 
describing flux within the rod, the bessel function Io(Wor), 
and flux within the moderator AIo(Wir)+ BKo(Wir) to match 

the experimental flux traverses for the 0.75 inch rods as 
shown by FigureA-3. Brookhaven National Laboratory found 

it necessary to multiply the XK of the fuel rod by an average 
value of approximately 1.5 and the X of the moderator by 
1.25 to match theoretical with experimental flux traverses. 
When these adjustments to ® are made, the calculated value is 
within four percent of the measured value of the disadvantage 
factor for 0.75 inch rods. 

From this limited comparison it appeared that the one 
group diffusion theory with zero neutron leakage from the unit 
cell would give sufficiently accurate answers at very low 
enrichments and at diameters below 0.6 inch. However, as 
enrichment and rod diameter increase, the disadvantage factors 
become increasingly less accurate. For the 0.75 inch case 


the previously mentioned arbitrary factors could be applied; 


69% 
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but until more experimental work is done with enrichments 


above 3% larger, undetermined errors will probably exist in 


the disadvantage factors computed by this method. 


Sample Calculation for Disadvantage Factor 
Stainless Steel Clad, UOo Fuel, Ho0 Moderator of 


~79 gm/ cm’ density 


ayy, ~ 

R, = Ro J Vijy42 

Vo = WP Ro* 

Vo = (.239)(R,#.019) 
mE = Ri? - (Vo + Ve) 
V, = Vuoo + Vo 


. 89 Vo/¥i iz. 
| Ezy, 1.60 
.212 (l-enrich.) .2065 
2.16 (enrich. ) 20432 
Zr. 2499 
ne i (3Ze 24, = 633 
x Hi Ro = .1675 
1/2 Wri = .290 
Fig. 3 (E-1) 006 
?, a 


Ro 


fe = 
KoRo 


Fig. 


= 1.036 


= 508 ecm 


0115 “H20/¥» 0106. 





21535 Vo/¥1 .OLL9_ 
Lo; 0225 
-O404(1l=-enrich. ) .0395_ 
9.65 (enrich.) 21932). 
Zz, 2326" 
13 2g, Zec = 415 


= .211 


1 F# 1.0054 
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Figure A-I 


Intracell Flux Distributions of 
Water-Uranium Lattices 
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TABLE A-I 
COMPARISON OF CALCULATED AND MEASURED 
DISADVANTAGE FACTORS 


0.6 | i308 | 3 [a | 1 | 5.7 | 
Tes | _usos [2 [aor | ase | 20 
Ce Ee Ee 
Tes | use | 9 [axe | use | 30 
a Ee 
Des [sass | as [oa [aver | or 
Toe | ascere | 5 [aes [axe [20 
0.75 | asco | 2 [2.93 | 2.595 | 26.6 
0.75 | asco | 2 [ase] 395 | a0 
ors | noose =| kt et | zp | eo 



















1. Kouts defines Vw/Vy as the ratio of water to uranium for 
0.6 inch rods and as the ratio of water plus clad to 
uranium for 0.75 inch rods. 


2. Using adjusted kappa values of 1.5 X for fuel and 1.25 4, 
for moderator. 
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TABLE A-II 


0.6 Inch Diameter, 1.15% Enrichment, Aluminum Clad Uranium 
Rods in Light Water (BNL 1987) 


oe) 2 | 6] | le 
Pe | ager | astoy | a.eso | -000 
= ase | aaen | a.zse | 050 


~~ 













lL. Dividing the moderator into cladding 


=, = .4h5 and coolant, thermal utilization is 
defined as 

Se, - .0195 - p, i / 

Eu = .0115 Pu Ey V, + Dy Zy V,, + Oy Z¢ Vat 





TABLE A-ITII 


0.6 Inch Diameter, 1.027% Enrichment, Aluminum Clad Uranium 
Rods in Light Water (BNL 1796) 
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ee 
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0195 
Eu = .0115 
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TABLE A-IV 


0.6 Inch Diameter, 1.3% Enriched, Aluminum Clad Uranium 
Rods in Light Water (BNL 1783) Values from graph. 
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APPENDIX B 
METHOD OF ADJUSTING HOMOGENEOUS PROGRAM CONSTANTS TO OBTAIN 
CRITICAL MASS OF A HETEROGENEOUS REACTOR 


The disadvantage factor was utilized to adjust homogen- 
eous program constants to give a critical mass of a heterogen- 
eous reactor. Each value for a thermal cross section is ad- 
justed to have the same value as the effective volume average 
value of the corresponding cross section in the heterogeneous 
reactor weighted by the relative flux throughout the unit cell 


in the following manner. 


Z,(25) = aera ee Vo Yo = N(25)@-4(25) 
Vi ) : 


where N(25) are the atoms per cubic centimeter of fuel and 


N'(25) are the atoms per cubic centimeter of unit cell volume. 


N4(25) = N(25) No 
Vig + ¥5. 
Therefore ' 
Vv Vv 
/ O 1 
2 Se 2 
o6(25) (*2 n =z) on ' 5) 
0 
This factor was used to adjust both U-235 and U-238 absorption 
and transport cross sections. 


Light water cross sections were adjusted in an analogous 


manner a , 
S,(w) = N(wi7a (w) Vz = N'(w)Ga (w) 
YvoPotVid 
N'(w) = Nw) Vz 
Vo Jae 
+ *¥ 
Ga (w) = / 2 oa w) 
(#2 - + z) 
i 


This factor was used to adjust water absorption and transport 
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cross sections. 

Cladding cross sections were not adjusted because the 
flux in the clad is very nearly equal to the average flux 
in the unit cell. This is apparent by inspection of Figure 
A-I. This i ys was checked by calculation by i 
values of W/ Save with experimental values of "Fu 
where 

&, = Vo #Vy 


Ga Vo % + Vi 
a 


The error in all cases indicates that by not adjusting the 


eross sections for cladding a greater cladding absorption 
was being assumed. This was a result of the fact that the 
average flux in a unit cell was higher than the actual 
neutron flux in the clad. However, the error introduced is 
necessarily small because of the small percentage of volume 


that the clad occupies in a unit cell. 


vl. 
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APPENDIX ¢ 
ADJUSTMENTS TO NUCLEAR CONSTANTS 
FOR POWELL FOURIER PROGRAM 

The Powell program was written for a highly enriched, 
homogeneous, heavy water moderated reactor. The adjustment 
of the cross sections to obtain a heterogeneous critical mass 
is discussed in Appendix B. Table E-I lists the inputs on 
the data tape for Powell's cylindrical Fourier progran. 

In this thesis it was planned to utilize the program 
for low enrichment, light water power reactors. However, there 
was no provision for the handling of uranium-238 in the pro- 
gram. The best method appeared to be to continue to define 
7 as J aa and to treat the uranium-238 as an addi- 

a 


tional cladding. By this method the absorption in the 


uranium=238 would appear in the four factor formula in the 


 # 


thermal utilization. 


epee 
the rmal ee in U-235 
total thermal absorption in core 


The fast fission factor, £, and the resonance escape proba~ 


f 


bility, p, were combined with the 7 for uranium-235. This 
product Expy with an additional arbitrary factor of 1.1 to 
bring experimental and program critical masses into agreement, 
constituted the input for "Pas" in the Powell program. 

The program, as written, included density and atomic 
weight for aluminum and heavy water as the cladding and modera-~ 


tor, respectively. In order to use the program for other materials, 
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it was necessary to either modify inputs to account for these 

constants or to modify the program by changing the inputs of 

the four drum addresses concerned. Both methods were utilized. 
Method A combined the effect of uranium-238 and clad 


material in the following manner 











(3) 









or (@?) = 228 Ve | Boy ¥ a 

VotVe) (6-025 x 10**) 

av ” \ “ : 

oe (at) = Za(¢) A_ =. ‘SZ. fab) 
4 (2) 6.0aSxl0"7 6,.0a5% 10"? = 


2 (20) = | er) ox (50) (5) 


The drum addresses were modified to read 
DA 420/ 420/ Pao [a0 
DA 42h/ 42k/ Aigo = + 2 
DA 418/ 418/ fat = +1. 
DA 422/ 422/ Aas 


ob das 


This method used the macroscopic cross sections divided 
by Avogadro's number and modified the density and the atomic 
number of the elad to equal one so that when the computed 
N(al) and "o~(al)" were multiplied together an equivalent 
clad macrossopic cross section was obtained. Transport cross 
sections were computed identically. 

Method B is a mathematically equivalent method and was 
devised to circumvent the modifying of drum addresses. this 


was done in the following manner: 





a 
—Elaee 
=> 
= 
é 
& 
=> § 
oe 
= e 
/- 
>_> 
> 
=> > 
= 
- = &S SS 
ja 
_ — 
= 
= 
a> 
— 
ad oe - 
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Vt V, 
‘Oz (Ql) = ~ Cac + Ox, eV, % | (6) 
Nat (Ve +l) 
(rg <1) 
" ‘ a 
ox (0,0) = he ! +h, Ce (H, 0) (7) 
520 
23 23 
ae = aes 4 2 110 |, JOS (3) 


Method B is merely a variation of Method A allowing the 
use of the program without modification to existing constants 
incorporated in the writing of the program. Method A, when 
used in the program, gave inconsistent low critical masses 
which led the authors to believe that the drum addresses were 
net being properly modified. This was confirmed by the sub- 
sidiary outputs of the program. The following page is a sample 
ealculation for UOo fuel and light water moderator of density 
0.79 gm/em3 employing method B. 

The quantity “x" in the program was changed from the 
ratio of the velumes of aluminum to heavy water to the ratio 
of the sum of the uranium=-238 and cladding material to light 
water. "x" is used in the program to homogenize the macroscopic 


cross section of cladding and moderator. 


z 


ay 


£ (w) at “8 (a) 2 + (9) 


NA 
rT 


Zt] + Zal}f eat (10) 








Sample Calculations for 
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Nuclear Data UO2 Fuel, H20 Moderator, Pp. = .79 (Method B) 


Lia Zin 
28 ~ O404 S22 
55 « 535 890 
Vo = .810 
V, = , 126 
V+ Vy = 936 


Vo + Ve x-06025 = .0564 (A) 


y) = 
= = 1.036 V5 xi = 
9, 
0 
Vot Vi_ 2.433 2 
Vo ~ vi, Zo49 | 
. % 

Vo + —B s 2h | 
ee, eat a 
Vote + Vo Zo 0 

4, 


-O404 x Vo x (B) ee, 








0 1535 x Vo = ED ad wit 

(A) = (0564 mf 
0212 x Vo x (B) =. 
©890 x Ve = .112 

» 280 il 
(A) = .0564 
344 x (C) = .347 ~ 
50.0 x (CG) = 50.5 = 
x 3 (« 8 al + 2126 = ° 21 
° g lr; 9 


e910 


497 
"oY Do0" 
"CZ D20" 


=, 616 


bf 
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APPENDIX D 


In order to account for the fissions in uranium=-235 
and uranium=-238 produced by fast neutrons before slowing 
down to thermal energies, a fast fission factor is utilized. 
It is defined as the ratio of the total number of fast nev- 
trons produced by fissions due to neutrons of all energies 
to the number resulting from thermal neutron fissions. Com- 
puted fast fission factors were calculated by using Hellen's 
empirical formulation: 


€ =1F idle 





+.288 JS (1) 


density of water 


a: 
<. 


Vw = volume of water 


Vy = volume of uranium 


iG 


Veg = volume of structural materials 

This factor was checked for eight critical experimental 
eases performed at the Bettis plant of the Westinghouse Cerp. 
for which a summary is given in Fable D-I. Errors of computed 
values are in all cases less than 1 percent of measured values. 

That fraction of source neutrons whieh escape capture 
while being slewed down to thermal energy is represented by 
the resonance eseape probability, p. The method outlined in 


Glasstone and Edlund for heterogeneous reactors was used with 


the empirical quantities from The Reactor Handbook, Volume I, 
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to compute resonance eseape probability. 

This method assumes that neutrons within the energy range 
where resonance absorption by fuel is important can be repre- 
sented as a single energy group, using effective slowing down 
and absorption eross sections. 

A resonance utilization analogous to the thermal utiliza- 


tion is defined as: 


". = rate of absorption of resonance neutrons (2) 


rate of production of resonance neutrons 


~~ 2 Ye (Qle Zao 


Equatien (3) assumes that the leakage in the resonance group 


(3) 


is small compared to abserption and slowing down. Resonance 


utilization may be expressed in terms of lattice parameters 





as: 
1 ig MZ F + (E-1) (4) 
r Vo 
ao 
where F and E for cylindrical geometry are defined in Appendix 
A. 
Resonance escape probability is written in terms of 
fy, as: 
Ir 
= e <> 
P ei ee (5) 


It has been assumed that the resonance neutrons can be 
approximated by one energy group. They exist over an energy 
range greater than would be predicted by one group diffusion 


theory because of the selective absorption of neutrons as they 
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pass through uranium=-238 resonances enroute to the interior 
of the fuel rod. 

To determine this effect average slowing down and ab- 
sorption cross sections are computed based upon empirical 
relations derived from experiment. The microscopic average 
resonance cross section is expressed as: 


stu § 


ord 
In E2/E, 


(#] = 


= mass of fuel lump, grams 


= surface area of lump, om* 


M 
s 
A = empirical constant, barns 
M = empirical constant, gm/em* 


Ei = lower energy limit of resonance 


E> = upper energy limit of resonance 


Table D-II lists empirical constants from The Reactor Hand-= 
book, Volume I, for natural uraniun. 


For cylindrical rods: 


Z,,- ma [it a He] - NA [i+ Sons 
a in” ~«Eo7E} n b2/5y 


The effective resonance slowing down cross section for 


the moderator may be expressed by: 


Zz, 7 af Fo 


Table D-III lists the resonance constants for moderator from 
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The Reactor Handbeok, Volume I. 

A sample calculation sheet for the calculation of reson- 
ance escape probability for UOo fuel is shown on the follow- 
ing page. 

Eight experimental cases wae ealculated and comparison 
of measured and computed values are shown in Table D-I. The 


maximum percent error from the measured p is 2.44. 


SS & = @& @& > 
= a, 
7 = a * - 
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Sample Calculation for 
Resonance Escape Probability UO Fuel, Ho0 Moderator (.79 gm/cw ) 
Core No. 39 


Fuel Rod Radius, R, = .508 em 





Unit Cell Radius,R, = -880 cm 
"282 = .753 
Ro 
+ 1.0 = "7 
“Toyez * «+9353 = +0620 Zao 
42 x Ro = .213 (AR) 


~649x Ro = .329 (W, Ro) 
.649xR, = -570 (%, F) 


From Appendix A, F = 1.005 
E = 1.026 


—“x ‘=—xF = AAO KI ge Oe 4,22) (A) 
Vo 2 os e 0 x ~O 


x 
fr= [tat sl = fh.22 + 1.026) “2 = .190 


f = ~ieoo - 
cr. ~ ‘eto * 238 (8) 
=mé = .790 


Fast Fission Factor 


ie » 1565 7 


W 
+875 Py HPO 
2156 = 1.041 


=1+ hoy .288) (.126 
1 + (.875)(.79) Ai ; 10) * Rey Ce 75)(. 1) 


L. Dap (1.1)(1.041)(.790)(2.08) = 1.885 


8 9) 
Ut 
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TABLE D-I 
COMPARISON OF MEASURED AND CALCULATED FAST FISSION 


FACTORS AND RESONANCE ESCAPE PROBABILITIES 

















1.043 | 1.047 
1.056 
1.066 
1.043 
1.056 
1.041 
5/1 1.029 


4/1 





TABLE D-II 
- 
EMPIRICAL RESONANCE CONSTANTS FOR FUEL 


[Material | A(barns) |/y (gm/om2) | in 52/81 












1. Reactor Handbook, Vel. TI 


TABLE D-III 
EMPIRIGAL RESONANCE CONSTANTS FOR MODERATOR — 


Material cet S&H 
(gm/cme ) barns/atom 





- Reactor Handbook, Vol. I. 


: Ae is multiplied by .88 for UO, fuel. 


Le 
2: 
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APPENDIX E 
POWELL CYLINDRICAL FOURIER PROGRAM INPUTS AND OUTPUTS 





Figure E-I defines the dimensions for the six radial 
and five longitudinal regions for Powell's program. [In this 
thesis two radial regions; a central core and radial reflec- 
tor and three longitudinal regions; a central core and a 
top and bottom reflector were used. Thus R, and R2 were 
effectively set equal to zero, i.e. 1 xX 10729 em, R3 was the 
core radius computed from heat transfer considerations. The 


thiekness of the two outer reflectors was set equal to zero 


by letting R3 = Ry = R5. Re, the outer radius of the reactor, 
was then simply the sum of the core radius and reflector 
thiekness. In every case it was assumed that a reflector 
thickness of eighteen centimeters was equivalent to an in- 
finite reflector. The radial reflector thickness was set 
equal to eighteen centimeters. The thickness of the top and 
bottom reflectors (number two) was set equal to zere. The top 
and bottom reflector (number one) thickness was set equal to 
eighteen centimeters. 

For the radial region chosen as the eore, in this ease 
region three; Fermi age, thermal diffusion coefficient, in- 
finite multiplication factor, and thermal absorption cross 
sections are computed. Thus, the inputs to these regions 
are normally chosen as zero except for the thermal diffusion 
coefficient which appears in the denominator and must not have a 
zero value. This same situation exists for these quantities 


in the longitudinal properties. The core properties are 








= 
==> 
Sh <n 
in 
7 oe = 
—> pa 
SS 
_ = 
| 
<-> _ 
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replaced by computed numbers as the program runs so that any 
number used to deseribe the core region will be erased. 

The inputs in the nuclear constant section are discussed 
in Appendix C. 

The final inputs are two kop guesses which are used to 
compute values of the critical equation. In general, neither 
of these guesses will satisfy the critical conditions, so 
the two assumed values of Koo correspending to the two com- 
puted values of the critical determinant are extrapolated 
linearly to the value of Koo at which the value of the deter- 
minant is zero. ‘his process is continued always throwing 
away the guess that is most in error until the difference be- — 
tween M(25)n =3 and M(25), =2 is less than one tenth of a gran. 
The guesses need not be close to the answer except that the 
closer they are, the less time will be required for convergeice. 

Of particular interest when making a parametrie study 
such as this is the fact that once the main program is fed into 
the computer, a series of data tapes may be run without read- 
ing in the main program again. The data tapes need not be 
complete, but need only be cut from the beginning through the 
numbers which were changed from the preceding data tape. All 
constants read in on the first parameter remain in storage 
until replaced by a new number, and all computed values in 
storage are erased prior to starting a new run. Thus when 
systematically varying one constant a great amount of time 
can be saved. 


Table E-I is a sample data tape preparation form from 












the J. Powell thesis using data from run 39. Table E-II is 
2 sample of the form found convenient for typing this data 
n the Whirlwind data tape. 
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Figure E-1I 
Powell Fourier Program 
Core Configuration 





le-—R, 








TABLE E-I 
DATA TAPE PREPARATION FORM 
FOURIER METHOD FOR CYLINDER 





fe TAPE 404-307-39 LEE AND GARDNER 
(24,6) 
NOT PA 
///Dimensions 
DA 32/32/ + 0. 
41.0 x 10719 Ry 
+1.0 x 10729 Ro 
+ 94.5 R3 (eore radius) 
+ 94.5 Ry 
+94. 5 Re 
+112.5 R6 (outer radius) 
+243. Core thickness 
+18. Thickness top reflector 1 
+0. Thickness tep reflector 2 
+18. Thickness bottom reflector 1 
+ 0. Thickness bottom reflector 2 





/// Radial properties 


/// Thermal absorption cross sections 


+0. Zal 
+0. Bar 
+0. Z a3 


0. = ah 
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| 
‘ 
°*ea @F & t 


(oe oom, ° 
a 
> eel? Oe Oa PS 





> iy ee” ee ee eee 


— 6 ee 


'_ -_ es ese ao | —_—€ mie 








+ 0. Fad 
+_.01145 E ab 


//1 Xoo + 0. Led 
+0. ky 2 
+0. ko 3 
+6. ko 4 
+0. ke 5 
+0. Kyo 6 


/// Thermal diffusion coefficients (non zero values) 
+1. D) 
+l. Do 
+l. D3 
+ 1. Dy 


*~- 1. De 
$ 62 De 
/// Ages 
+ 0. x1 
+8, WY 2 
+ 0. Tq 3 
+0. 7 4 
+ 0. A 5 
$49.5 7 6 
/// Nuclear constants 
+ 421. Oa (25) (barns) 
+ 350. co” f (25) (barns) 


+ 1.885 % (25) 
+ 2.46 y (25) 
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+ .910 O~a (Al) (barns) 
f+ _. 347 ona (Do0) (barns) 
+ -616 xX = Vol. Al/ Vol. D0 
2:4) 3 4,2 %+ i hee 
+0. th 2 
epee 
+0. f 
150.5 7F-(0,0) 
set BY a Oy (Al) 


+1. Fermi fraction of age 
ro. Diffusion frac. of age 
t 0. 


/// Longitudinal Properties 


/11 Xoo 


+t .01145 Za Bottom refleetor 2 
+ .01145 Sa Bottom reflector 1 
+0. 

+ .01145 Sa Top reflector 1 

t .01145 =a Top reflector 2 


+0. ko b2 
+0. ko bl 
+ 0. 

+ 0. aS #8 





+ QO. k,o 2 


/// Thermal diffusion coefficients 


+ .2 D 41 









/// Ages 











Note: All numbers must have #¢ 


or 


- and decimal points. 
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TABLE E-II 
SAMPLE DATA TAPE PRINT OUT - RUN 39 
fe TAPE 404-306-39 LEE AND GARDNER 
(24,6) 
NOT PA 
DA 32/32/ +0. #1.0 x 10719 +1.0 x 10710 
$94.5 494.5 494.5 #112.5 +243. 
+18. +0. +18. +0. 
+O. +0. + 0. +0. +0. +.01145 
+ 0. +0. +0. +0. +0, +0, 
+1. +1. +1. en #1. + .2 
a0. +0. + 0. +0. +0. t49.5 
+421. +350. + 1.885 + 2.46 + .910 + .347 
+ 6616 449.5 +0. +0. +0. $50.5 
$4.97 4¢1.0 #0. 70: 
+ O1W5 +,01145 4.01145 4.01145 +,01145 
+0. +0. +0. +0. 0. 
Pe + .2 + .2 + .2 + .2 
th9,5 +49.5 +49.5 +49.5 
+1. +1. +1.d. $1.2 +0. +0. 


4 START AT 186 
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OUTPUTS OF POWELL CYLINDRICAL FOURIER PROGRAM 


The program used in this thesis computed critical mass 
and infinite multiplication factor. However, the printout 
format of, the program includes several other numbers which 
ean be very useful in the evaluation of the answer received 
and also very confusing to the uninitiated. 

An actual printout received for run 39 is shown below: 
TAPE 404-306-39 LEE AND GARDNER O541 2 4-16-57 
DECIMAL 

+.13461303/=00 +.17050170/~00 
+. 30481165/+06 +.30373265/*06 
~.15789275/=01 = .72544313/+02 
¢ -11931752/ O1 #.11916215/+01 


+ .30481165/ 06 -.15789275/=01 
+. 30373274/ 06 - .70792410/+03 
+.16155871/=27 +.16155871/=-27 +.16155871/-27 +.68174107/+07 
+.29004734/-00 +.49500001/+02 +.49500001/+02 +.49500001/+02 
+.16155871/=27 +.16155871/=27 
+.23188179/=-01 +.20502647/-01 


Interpretation would be as follows: 

In the tape number the 404 is the problem number, the 306 the 
programmer s number, and the 39 the run number. The numbers 
fellowing the name are the time and date of the run. In each 
number printed, the numbers after the slant (/) lecate the 

decimal point of the preceding number (+.11916215/#01 = 1.1916215). 
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The number +.16155871/-27 signifies zero. The first eight 


numbers are the following: 


nR=2 a= 3 
Time started Not used Not used 
Critical mass-M(25) (grams) 304811. 65 303732.65 
Value of critical determinant- A ~.015789275 -.0072544313 
Infinite multiplication factor-£, 1.1931752 1.1916215 


The next four numbers represent the approximation before 
the final number printed out above for the n = 2 and the n = 3 


terms in the truncated Fourier series. 


n= 2 (M(25)) es =em (A) 
i-1 approx. 304811.65 em. i-l approx. = ~.015789275 
n = 3 (M(25)’ n=3 (4) 
i-l approx. = 303732.74 em. i-l approx. = -.0007079241 


The following ten numbers are read across the page from 

left to right and signify the following: 

The first three always remain .16155871 x 107£7 = 0 

(4) core volume -= 6817410.7 cm? 

(5) core thermal diffusion coefficient -- .29004734 cm. 

(6) core age == 49.500001 cm@ 

(7) average reactor age -- 49.500001 om? 

(8) average Fermi age <-- 49.500001 om 

(9) average diffusion age -- .16174107 x 107¢7 = 0 

(10) This number always remains zero signified by .16174107 x 107¢7 
The last two terms are the number of atoms per cubie centi- 


24 


meter x 10. of clad and moderator, respeetively. 
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(1) clad -= .023188179 atoms/cm? x 10724. 
(2) moderator -- .020502647 atoms/cm x 10724 

With these outputs it is possible to compare eaeh of 
the outputs with inputs to determine whether the program is 
operating satisfactorily. In addition by comparison of the 
i-l and i approximations for each of the terms in the trun- 


cated Fourier series, it is possible to establish hew fast 


the program is converging. 
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APPENDIX F 
NUCLEAR DATA 


A survey of the literature was undertaken to determine 
best values of the nuelear data that would be required. In 
general the best sources of information were found te be in 
the various Breokhaven National Laboratery reports and in the 
Geneva papers. 

A. Cross sections 
All cross sections were obtained from BNL=325 except 


where noted. 
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TABLE F=I 


CROSS SECTIONS (2200 m/sec) 
by 


— -2 
we lL & oF 2 ae 
Nuclide (barns) (barns) (barns) (bafns) (bafns ) atoms/cm 


enmepcmpecmemeriisermees 3 nee ee Cm eee CEE RR rer GE nt eee? 


H . 330 . 293 38.0 0334 
Ho0 . 660 - 586 | 0334 
U-235 687 599. 10r. 580. 503. O47 
U=238 2.75 2.44 8.3» » O474 
Al > 230 204 1.4 .06025 
Fe 2.53 2.2h ag . 0622 
Ni 4, 6 4.07 ay a5 00727 
Cr 209 2.57 3.0 0149 
Zr 19 168 8.0 ~  ,0423 


1. O— signifies effective cross section-corrected for Maxwell 
Boltzmann distributien. 


2. Type Nr. 304 stainless steel was used. 





In pe 


Density of Zr 


TABLE F=II 
CROSS SECTIONS (508°) 














_ Nuclide NEO <a 22, on ral 
U-238 0474 1.80 .0839 9.45 bho 
§-238(U02) -. 0224 1.80 .040% 9.45 220 
U-235 0474 = «431. 20. 96.5 4.57 
Fe 0622 i Se 5c 5 Ns ye 3) 
Ni -00727 3501 0219 20.0 » L46 
Cr 0149 1.90 © 0283 3.34 ~050 
Stain S. £535 890 
Al -06025 elL51l .0091 1.365 0823 
[0 0264 432 .01145 62.9°> 1.66 
ar 6. 0423 1245 .00526 7.95 336 
- = Og Nr | 293-2 = 65502 

a 
2. 0g = (1G) (1-W/5-d; Ao= 3h 


Experimental=-Hughes, Pile Neutron Research, page 223. 





Density of water = 0.79 gm/em? (508°F) 
Density of UO02 


10.07 gm/em3 
6.4 gm/em 








B. Fermi Age 


Early measurements of Fermi Age were carried out at Oak 
Ridge in 1944-45, for water to metal ratios in the range from 
1 to 3.5, using rod diameters of 1.18, 1.10, and .787 inches 
and enrichments of the order of 1%. A fission neutron source 
was used, and the rods were wrapped in cadmium to prevent fur-= 
ther fissions in the lattice. The activities of cadmium cov= 
ered indium foils spaced through the lattice were measured, 
and age to indium resonance computed by 

v2 ke 


A correction term for the difference between indium resonance 


a 


and thermal was added, but this is small (about 1 om®) , 


These measurements gave a set of characteristic curves 


‘ showing an age of about 40 on” 


at a volume ratio of 3 ¥o 1, 
rising sharply as the volume ratio was reduced below 1.5.. Mm 
increase also occurred with reduction of the rod diameter. 
Kouts, in Geneva Paper No. 600, “Exponential Experiments 
in Light Water", has reassessed these data in terms of the pos-= 
sibility of distortion through fast and epithermal fissions. 
These would not be prevented by the cadmium wrappings and 
would tend to occur more strongly at the lower volume ratios. 
They would have the effect of increasing the apparent age. 
_ Kouts states that, in fact, the fraction of fissions produced 
by neutrons above cadmium resonances (0.4 ev) at a volume ratio 
of 3 is about .11 while at a volume ratio of about 1.5 it 
rises to .58: this includes both fast fission in U=-238 and 
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epithermal in U-235. 

Kouts made measurements of Me the migration area, in. 
low enrichment water lattices with volume ratios of 1 to 3. 

The age computations which he felt were the most accurate 
were based on the relationship 

A > eget = exp (BM) 
His method was that of poisoning the lattices with increasing 
amounts of boric acid, computing B* by differentiation of the 
thermal flux (found by measuring indium foil activities), com- 
puting f, and solving for M“, considering the other factors 
to remain constant. 

His results were markedly different from the Oak Ridge 
curves. They show a curve very slowly rising with decreasing 
volume ratios, and, indeed, in one case flat within the experi- 
mental error. Measurements made with 0.600 inch rods are con- 
sistentl greater than those with 0.750 inch rods, but can still 
be considered almost flat, with a mean value of about 30 om, 

Chernick discusses in his paper, read at Geneva (No. 603, 
"The Theory of Uranium Water Lattices") methods of calculating ~ 
age, and, while he does not show actual calculations or tabu- 


late results, he does give a simplified formula 


M2 

1 ou 

T = 3 [a + Bure ji 

which does not include volume ratio as a variable and states 
that “Calculated values of the age of fission neutrons to 
indium resonance in water to uranium mixtures show little 


variation if the uranium to water volume ratio is less than 
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unity...,For equal volumes of uranium and water, calculated 
ages have run about 20% higher than the value for pure water." 
Chernick gives as the best current value for age to 
indium resonance in pure water, 25 em* , 
1.2 x 25 = 30 cm 
Also included in Chernick's paper is a formula attri- 


buted to Ehrlich and Deutsch of KAPL 


y A 7 iin0 (1+x)e 


(P+vx) 
where v = volume ratio, uranium to water 


= density of water 


+ 


4 
8 


0.38 for stainless steel 
0.40 for Zr. 


FT = 26 (tay. 52.5 om? 


Deutsch in Nucleonics (January 1957) gives a method for 
computing age in mixtures of waters and various metals such 
as Al, Fe, and Zr, for high enrichment reactors. These ages 
will be in the range of 40 to 50 or more, but they are appar- 
ently based on significant volumes of cladding material and 
insignificant volumes of uranium. These, then, do not fit our 
problem. Because of the similarity between this method and 
the formula ascribed to Deutsch by Chernick it seems possible 
that that formula might well have been intended for use with 
high enrichments also, and the "v" in it might have been in- 
tended to express total metal to moderator ratio, being largely 


made up of cladding volume. 
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Finally, J. E. Wilkins et. al., in Geneva Paper No. 597, 
"Status of Experimental and Theoretical Information on Neutron 
Slowing Down Distributions in Hydrogenous Media", discuss com=- 
putation methods for strongly absorbing media in which thermal 
diffusion area is small compared to age, and compares results 
to Kouts' data. The calculated line follows Kouts! points 
for 0.600 inch rods closely but turns sharply up at a volume 
ratio of unity to a value of almost 38, compared to Kouts? 
measurement of 32.5. 

Since Powell's program offers a choice of using the dif- 
fusion or the Fermi formula for the age, for our purposes 
it appeared best to use the data of Kouts! and Chernick's theory; 
since Kouts! data was obtained from the physical reality by 
the Fermi formula, the machine computation will also be based 
on the Fermi formula. 

The value selected was a constant value of 30 om? age to 
indium resonance for water to uranium ratios of 1 to 3. A 
1 em* correction to thermal from indium resonance gives an 


age of 31 em* for all enrichments and rod sizes. 
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C. Estimation of reflector savings 
Theoretically, savings for an infinite water refleetor, 


sheuld appreach 


| = » L 


D r 


3 


Sinee o¢ for water is se large, it will overshadow om for y 
and SS, and 


Bm = 8: 


At room temperature, L for water is 


= ke ow 


ee ee 
os 
However experimentally, values range from approximately six 
to eight. L will also vary with temperature. For the econdi- 
tion in this thesis 


thes - | Ze Zz = Neo | Ox cy 
L cold Fi, Sy, H Hy CEH 


1 
a ple == i 

5g X 1-55 ? 

We assume that the refleetor savings will vary with tempera- 


ture in the same way that L dees 


An = 1.7 

A 
By seleeting the end of the experimental range closest te the 
theoretical prediction 


dx = ar vg 21:7 x6 =10-2 & 10 en. 


a —_ . a a SS 
_ 
e = FF 

> > =z: 

se ' 
ee 

7 

a = => 

- @ '. 
= ® 
* 
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APPENDIX G 
COMPUTATION OF CORE RADIUS AND NUMBER OF RODS 


A methed for determining necessary core dimensions fer 
the transfer of predetermined thermal powers, given the unit 
cell dimensiens, was developed from Rehsenow, Lewins, and 
Barger's paper, "Steady state Temperature Bia entees in a 
Nuclear Reactor with End and Center Fed Ceolant."  Rehdewen, 
Lewins, and Barger develep functions P and Q for the tempera- 
tures of the fuel centerline, the fuel clad interface, the 
eceolant clad interface, and the eeolant bulk temperature. 
Function P expresses, in dimensienless form, the temperature 
differenee between the point in question and the ceolant in- 
let; function Q expresses, also in dimensionless form, the 
heat transfer resistance between the area and the eoolant. 

It is shown in the paper that the maximum P along a chan- 
nel, and hence the highest temperature, is 


P= Q?+l 

Since we were investigating maximum performance of 
power reactors, the limitations placed on our cores were 
principally those of temperature, and, of these, it appeared 
that temperature at the fuel rod centerline and of the coolaht 
bulk temperature would be the most critical. Further, it ap- 
peared that if each of these Limitations were just met at 
full power, that the most efficient reactor, under the pre- 
scribed conditions, had been designed. 


The functions given by Rohsenow, Lewins, and Barger for 


1. Rohsenow, Lewins, and Barger, op cit. 
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the end fed case are ag follows: 


Coolant Bulk Temperature: 


where At, = t coolant max = t coolant inlet, OF 


Gn = heat transfer rate at midpoint of rod, piu ee? 


Ww = weight rate of coolant flow, 1b./hr. 
c = coolant specific heat, Bae 


A = channel cross section area ft - 


ag = rod heat transfer area, midpoint to end of core, rt" 


ap = rod heat manganese” area extended through reflector 
savings, ft 


Centerline Temperature: 


P=m4't 5 We “0 - sin wm 20 
On A ar ar 

ge Tw. fo |t,2,Rey 2 x lo (%%o) - 1 
A ar Ue Ro Ake I, (HoRo) 


where h = film heat transfer coefficient 


Ue = effective heat transfer coefficient for bond, clad 
and scale 


Ro = radius of fuel region 
R, = radius to outside of cladding 
1  - Rg 1n(Re/Ro) 1 1 
a eee) a te 
the third term expresses the resistance in the rod, assuming 


that 
@ = ATo (KoRo) 


Ate=t of fuel centerline - t coolant inlet 
¥ 
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In order to obtain a simultaneous solution for given 


At, At, determined by fuel melting and coolant 
polling t's 


h, Ue 


0, Kp, kg, © determined by materials and dimensions 
of the unit cell 


3 


in each case, Pp was not equal to Q-+1, In the equation for 
the coolant bulk temperature thus obtained, this leads to 


P= 1, and the resulting equation was solved for w 


w= (L+ein P=2)(a ai) 
T Ate -c. 80/ar 
Physically, this means that we have set coolant flow suffici- 
ent to maintain the bulk temperature within the limits set. 
This value of W was then substituted in the fuel center- 
line equation; that is, knowing the coolant flow rate, we de- 
termined the amount of heat that could be generated in the rod 


without exceeding the temperature limit set. 


ao 2 2 
r Atr* . (1+sin g ae) (Aan) - a. (32) 
qth a7? 
T At. e~(8q/ar) 


a 

oO 
/ ene Sy , + sin ¥ ap) (Aqn ) c %% | 
qn T O*e- c. 8o/ar A af a af 


a, 2 | 2 
= 120% (1¢sin J ge) (agh)?. (22) ar + + 
2 tel 


TT Dif - om (80/az ) “tf 
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a a 2 2 
te (l+sin F 0) in WT 2. ~ (At 
qt? = L : ar be AM, Oe Zz ar | (At,) 


a 2 2 
(1+ sin — Xx Qp 


2 
where one [ied + o. lL , PWo(H¥ Ro) -1 
o Ro” ke )=6 THR) 

For a particular series of reactors in which the rela- 
tion between core height and reflector savings remains the 
same and in which the same temperature limitations apply, 
this can be expressed as 

ager te / Sq 


2 a7 2 2 
were Tt - Jate(atein F ne At, sin 7 20 | (4*<) 
a 2 
(itsin J 2) 


In the class of reactors under consideration, core heights 





will be of the order of 8 feet and reflector savings of the 


order of 10 cn. 


a 
in (Ir coe = gi - da = “O- a, 
S - 8 n(L 113) sin 83.0 0.991 


9 


of sin q ar is taken as constant and equal to 1, 


T~ = Ate (Ate - Ata) 


q*% is the heat transfer rate at the midpoint of the 


central rod. The heat transferred in the entire rod, q,, is 
ic 
ao = at f cos (7 a 
~8o oi ar) da 
boa 


= ee ak cin (x =e = Me 
= 
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And the average power developed over all the rods, (assuming 


4 Jo flux distribution) q, is 


R 
Sree _ fa, J 2.405r. 2Trr ar 
re 
R! 


| R? 2. 405R 
q “90° 3H B- J (24 ) 


where: R €& core radius 


R +) 


wv 
at 


— 
f 


= reflector savings 


The thermal power generated in the reactor, 


P, = a xN where N = no. of rods 
Pa =q = 40x - 83 a Jy es R) 
N R 
— .837Td5 R(R+) ) 2.405R A,, = cross section area 
Pp te 45, oer of unit cell 
U 


This may be solved for R, and the number of rods may also 


be found from 


N= WR“/A,y 

Hot channel'factors are engineering safety factors applied 
at appropriate points in the heat fransfer calculations where 
uncertainties exist to assure a safe design. They may also 
be used to correct for the non-uniformity of neutron flux 
across the reactor, but this was accounted for in the process 
of averaging an assumed Bessel function Jo flux distribution 


across the core to obtain the average power developed in one 
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rod. We used only the engineering hot channel factors, there- 
fore, and adopted a set very similar to those described by 
J. W. Simpson, et. al., in "Description of the Pressurized 
Water Reactor Power Plant at Shippingpert, Pa.", Geneva Paper 
No. 815, which were as follows: 

ls > an allowance for the uncertainty in the heat. 


transfer coefficient, h, applied to reduce h. 


an allowance for the uncertainty in the allowable 


Cx) 


FQ 
heat generation rate in the rods; applied to re~ 


duce qn . 

Fa ry 3 an allowance to assure that the temperature 
rise stipulated is sufficient to prevent boil= 
ing, applied to reduce At,- 

The values of these factors are based on specified tol-= 

erances, known ranges of uncertainty, and engineering experi-~ 
ence and judgment. The values which we adopted for our heat 


transfer calculations were as follows: 


Fo = 1.53 
= 1s 
F At= 1.39 


By analysis similar to that shown above, the thermal 
fluxes required to produce burnout and local boiling can be 


computed and compared with qi 


2 a t 
Inbo " Aty) (A ty dno 
Burnout ¢ “T.2 oi we (Eel? 
h i, 
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it 2 Ate 
Boiling: AmB . Atwo (Atwo ~ F 


Values for thermal conductivities and the film coef-=- 


ficient were selected as follows: 


Btu oF 
Material k pr ft. 
Uranium oxide, U0. 0.95 
Stainless steel, type 304 10.6 
Zircalloy~=2 Sal5 
film coefficient h = 6000 BEU/L 64.2 OF 


A sample calculation of core radius follows. 
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Sample Calculations for 


Core Radius UO, Fuel, Stainless Steel Clad, Core No. 39 





tm = ty (Ate - Aty) = 4313 (4313-101) = 18.2 x 106 
7 = 4270 
2 O 

1 =— “ “4 ft pr F = 2.55 x 1077 
* nae 2.55 x 16 —_ 

R i 
me = Re 1” °/Ro . (.215) la (:336) = 1.22 x 107" 
uv Ko 12 x 10.6 
i), = .211 
Ro =. = Lo (2% Ro) =1 
Kr KR, Tv (&Ro) 

S22 eee ee _y 
1ox-5 sir * — 94.5 x 10 
Qe = 98.3 x 1074 
. a aa 4270 = 4Q. 10° 
‘mm = Qn Tex 98.3% 10-4 ms 
Burnout Check 
hed = 2.55 x10 + 1.22x107 = 3.77 x 10-4 
. 2060  _ 2060 = 547 x10? = 2 a" 
fe Bi = f5 y ° “Sm” . 
"BO ht f 


Local Boiling Check 


dn ee. 6 
Q 1.39 
aw 40.3 x 107 


1 § 


> 





ba 


1/2 core height, H = 4 ft. 


Radius to autside of clad, Ry = «215 In. 


Radius of unit cell, R, = .880 cm. 


. | 
do = af x 8R,H = 8(40. 3x10 pi-24 ee) 2.93 x 1077 


= .0676 mw 
2 
R(R +10) x Jy( $2058) = Te = 5510 


Core radius, R = 94.5 om. 


2 
Number of rods, N =. — = 11540 
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APPENDIX H 
CORRELATION OF THE POWELL FOURIER PROGRAM 
WITH EXPERIMENTAL DATA 





A series of test runs was made with the Powell Fourier 
program for Whirlwind, using data from experiments, in order 
to compare the results from the program with experimental 
results. The experimental data consisted of the six lattices 
tested by Krasik and Radkéweky— (aluntaee clad uranium metal 
fuel in light water, with characteristics as shown in Table 
H-I) and the three lattices tested by J. R. Brown et al. 
(aluminum clad uranium oxide, U05, lattices also in light. 
water), whose characteristics are also shown in Table H-I. 

Since values were published for the fast fission factor 
and the resonance escape probability for these lattices, 
these values were used in the computation of the program 
input "TEp". However, values were also calculated for ¢ 
and p by the methods shown in Appendix D. The largest error 
in the calculated "”E p" was 1.9%, and the average about 1%. 
Adjusted cross sections were computed as discussed in Appen=- 
dix © and in the Procedure. At first these were adjusted 
to an arbitrary density and atomic weight of one, and this 
value was sent into the computer to replace the values for 
aluminum which were a part of the program. Later it was 


found that these corrections were not being made by the 


1. Krasik and Radkowsky, op cit. 
2. Brown, Nooroff, Frantz, Volpe, and Harris, op. cit. 
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machine, and the cross sections were then adjusted to the 
aluminum values, which were left undisturbed. 

These runs gave results for critical mass of U=235 which 
were consistently too high by about 30%. In an effort to find 
a parameter which could be used to adjust the results to bring 
them closer to the experimental results, runs were made in 
which the input "NE p" was varied by a series of arbitrary 
factors. This is the input in which the greatest uncertainty 
exists; it is also one which may be changed without requiring 
any compensating changes in other inputs. And finally, the 
calculations are quite sensitive to it; a change of 1% in 
7EP changes the critical mass by about 2.5%. The results 
of this series are shown in Table H-II and Figure H-1 in which 
fractional variation of YEP from the calculated value is. 
plotted against the fractional variation of the resulting 
computed critical mass from the experimental. 

It will be seen that these results form a reasonable 
family of curves and appear to be fairly consistent. In order 
to give an average value of one for Mos (calculated)/Mgs (experi- 
mental), an arbitrary factor of 1.1 was chosen for Zep/VEp 
(calculated). This nethod of correction was chosen over the. 
alternate method of using 7é p(ealculated) and correcting the 
resulting Mes by a factor of 1/1.3 on the following basis: 

1. The curves were steeper at Vép/VHEp cale. = 1.1 
than at 7EP/ PEP calc. = 1.0. Thus an error of 
1% in pép would make an Mg5 error of only about 
2%, against an error of 5 or 6% if ep were in 





115. 


the other range. 

2. The spread of Mgcscalc./Mocexp. was narrower, 
about * 8%,at VE P/PE p(cale.) = 1.1 than that 
at VEP/YEP calc. = 1.0, where it was about 
* 10%. 

It thus appeared that both the error and the range of 
error would be smaller using the arbitrary factor of 1.1 in 
yep than if an arbitrary 1/1.3 were applied to critical 
masses computed frompep calc.. 

There is some measure of physical justification for this 
correction. In relatively tight lattices such as these, 
there is a considerable amount of epithermal fission, amount-= 
ing to as much as 10% or more of the total fissions, accord- 
ing to Kouts t+ If these did occur to a significant degree, 
they would tend to reduce the mass of U-235 necessary for | 
criticality, but would not be reflected in the calculation. 
And if they do occur, a not unreasonable way to account for 
them seemed to be by a fifth factor in Kyo » Which 1s effec- 


tively the correction which was made. 


1. Kouts, Price, Downes, Sher, and Walsh, op. cit. 
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TABLE H-I 
Results of Critical Experiments 


(%) Rod Nr. (25) 
V1/Vo Enr. | Diam. | Rods | x, 38 





5 I 


2 
1 
3 
2 
3 
3 
ny 
5 





1. Cores A, B, ©, D, B, and F are uranium fuel measurements 
from Krasik and Radkowsky. 


2. Cores 12, 13, 14 are UOo fuel from critical measurements 
by Brown, Noordoff, Frantz, Volpe and Harris and f, p, § 
measurements by Krantz, Smith, Klein, and Baer. 
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TABLE H-II 
Computed Results for Varied HEP 


M85) qe ¢ - 
(gm) x 10 x 104 







119. 


APPENDIX J 
COMPUTATION OF THE ADDITIONAL ENRICHMENT 
REQUIRED TO ATTAIN SPECIFIED CORE LIFE 





In the body of this thesis, a series of cores are con- 
pared which are designed to have equal power outputs; that 
is, they are eapable of the same rate of heat generation. 

In order to select a true optimum, however, the designs 
considered must be subject to a second requirement: they 
must be capable of the same total amount of heat production, 
or “power-time" of operation. The limitation on power-time 
is that the reactor must remain critical throughout the 
period of operation. 

The principal nuclear reactions affecting the reactivity 
of a low enrichment uranium fueled reactor are the following: 
U~235 -=- (n,f) FP 

a... (n,J) U-236 
U-236 -- (n, 2 (3 ) Pu-239 -- (n,f) FP 

: ----~ (ag’) Pe-2eo) —- Ge, Y) Pu-21 -- (nt) FP 

Production of the fissionable isotopes Pu-239 and Pu-241 tends 
to increase the reactivity, and their destruction to decrease 
1t as does destruction of U-235. Production of the neutron 
absorbing isotopes U-236 and Pu-240 and the neutron absorbing 
fission products (FP in the above expressions) tends to decrease 
reactivity. Since these production rates will be different 
in every reactor, depending on the quantities of U-235 and 


U-238 present, the flux at which the reactor is operated, and 


the resonance absorption probability in U-238, each just 
critical reactor will have a different, and relatively short, 


allowed time of operation before ultimate loss of criticality. 


In general, it will be necessary to increase the initial en- 
richment to a value greater than that required for just 
reaching initial criticality in order to attain practical 
power-times. This appendix develops a method for finding the 
factor by which the minimum enrichment for initial critical- 
ity of a given reactor must be increased to attain specified 
power-time. 

While the core life is expressed in power-time, it is 
useful also to define flux time as a measure of the degree 
of exposure of the fuel. The number of reactions which take 
place in a medium of given as which is exposed to a flux @ 
for a time T is 


Zp? 


Since, in the practical case, flux will be varied over a 
period of time due to changes in power demand, the total num- 


ber of reactions over such a period will be 


& |. 4 (T) ar 


The quantity [7 6 (T) af is defined as the flux time, t. 
© 


Following the method of Spinrad, Carter, and Egeler, 2+ 
the following nomenclature is defined: 


1. Spinrad, B. I., Carter, J. ©., and Eggler, C., Reactivity 
Changes and Reactivity Lifetimes of Fixed Fuel Elements 


120 . 


in Thermal Reactors", paper no. 835 presented at the Geneva 


Conference on the Peaceful Uses of Atomic Energy, 1955. 
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F, 28 = N28 028 eff where subscripts 25, 28, 
49, and FP represent U-235, 

as GO U-238, Pu-239, and pairs 
S25 N25 25 of fission products respec- 


tively. N, the number of 
Nig 9 n9 atoms per unit volume of 

the reactor, may be a func- 
tion of flux time. 


Z,49 

S pp = “pp Opp 
ICR, the initial conversion ratio, is the ratio of production 
of new fissionable atoms to the destruction of U-235 at 
startup: 


“an = Nog O 28 eff = 5.28 +442 25 25 € Py (1=p) 
25 ° 25 S.25 2.25 


where PS is the non-leakage probability during slowing down 
to the resonance energy. Since the effects of U-236, Pu-240, 
and Pu-241 are small they are neglected for simplicity. 

For the low enrichment case, differential equations may 


be set up for the flux time changes in the macroscopic cross 


sections: 
g 28 7 28 
d 2 oy S 
at 8,25 25 25 
so Sus = Tyo (pg - 2149) 
d S pp = ° pp Aes, _ S49 
at Lt Ags, 1 #Ang 


Where X is the ratio of capture to fission cross sections 
and the superscript "o* indicates the value at startup. Solu- 


tions of these equations are obtained as follows: 
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Ooct 
S26 = 126 —ala (2) 
-T yot 
E49 = S28 (l- ) (3) 


= ey 
Z, FP = A125 FP 1-e7° 25° t 4.28 |p ~FP (1-67 agt)|t 


1 * “os oO 1 +% ho hg 
Since Z 
fer = 0 
S28 aa = ICR x Jas (5) 


And equation (4) may be written as 


S, 26 14+%X 96 as 1 tH y 





Equation (7) expresses the total number of fission product 
pairs produced per atom of U-235 initially charged, at any 
flux time, t, and thus the total number of fissions, which 
is, in turn, directly proportional to the total heat produced. 

The productivity, // , 1s defined as the difference in 
rates of neutron production and loss per unit of flux time. 
When Tl = 0, the system is just critical. 

For each set of core parameters studied, the concentra- 
tions of all elements represented in the just critical, or 
“reference", condition may be found from the cell character- 


istics and the critical mass of U-235 as determined by the 





machine computation. For these reference systems (denoted 


by the subscript "R"), productivity is zero. 
T-0 = 07 = (25-1) 2 25_ -9 28_ ~51R ~ Leakagen = 0 


After the required flux time of operation, i] must again be 


zero. 
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(8) 


iT = 25-2) ei § lt) -ng-2) Aug - 228 7 2, bemeeae- GrP=0 (9) 


If the enrichment of the reference reactor is increased 


by a factor X 


aT = 257+) & S, 25p 5 J 28p - Sip - Leakegen 


Subtracting equation (8) from equation (10) 


Assuming that S28. $1, and Leakage do not change appreci- 
ably during the period of operation, equation (10) may be 


subtracted from equation (9): 
“a5-2) ( E\ 25-X |.25p) ~ (49-2) cing 9 5 pp = - OT 


2 FP as expressed in equation (6) may be broken up inte two 


parts 
‘el t 
SFP, = 2.25 = fice *® ) 
O omn t 
fms AF Se (oe “) 
lt Ang Ong 


Substituting values from equations (2), (3), (5), (13), and 


(10) 


(aay) 


(12) 


(13) 


(14) 
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(14) into equation (12): 


G2s5t t 
(M2571) (x S255 0 *? =x 2)25R)- XS 255 Opp (y__-° 25 ) 
1 tX%o6 O26 


~ 
$ ()ag-h) x A 25— ICR(1-e7 ™) 


a Snot 
Shot 4g 
_ 4g Ge Fe 
- ICR i" g-1) (1-e ~_ pet-s ) 
49 


= (25-1) A a5p(X-1) (15) 
Rearranging the terms and drepping out the common factor ee 


Oot 
weiss )(1-Noe- _ 1 Trp \4 ton }- CFPt 
?° TFxXos S25 Lug 





Sygt 1 Opp ) | 
+(1-e )(Yag-2 + TF X no Fig ? | 


=-X(125-1) t ("] 25-1) (16) 
Equation (16) gives the required factor of additional enrich- 
ment as a function of flux time and Initial Conversion Ratio. 
Appropriate values of the constants in equations (7) 


and (16) for the conditions specified for the studies are: 


O pp = 55 barns 


O26 = 437 x 25° 204 i 25 = 2.08 


Ong = 687 XK yo = -420 ug = 2-03 


Substituting these values, numerical forms of equations 
(7) and (16) are obtained: 


X = 


Np 
Noe 


125 


1.08 (17) 


1.08-1.19(1-e7 4324 roR [ 1.09(1-e7* ©87F)-. 0387¢ | 


_ (a-e7""97*)4 tor | .304¢-.4H1 (1-27 687t) | (18) 


Where t is flux time x 10741 


0 


Using these equations, X and —t may be plotted against 


25 


flux time for a series of ICR's, as in Figure l. 


fo) ss 
Since both ICR and Npp/N25 depend on X, in order to 


determine X from Figure J-1 an iterative process is necessary. 


The steps are: 


a le 


For the specified total number of fissions and 
from the critical mass of U-235 for the reference 
reactor, compute Npp/N28. 

From the unit cell geometry, enrichment, fast 
fission factor, resonance escape probability, 
and slowing down non-leakage, compute ICR as 

in equation (1). 

Enter the upper curves of Figure J=-1 with: 
Nrp/Nocp and ICRR, and find the flux time of 
operation needed. 

Using this flux time and ICRp, find X from the 
lower set of curves. This is the first approxi- 
mation. 


In the second approximation, enter the upper 


Fever J-1 
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set of curves with 


NFP/No? = 1/X(Npp/N2sR) 


a torn 4. = 1/%28\ 28 € P (1-p) 


These will give a different flux time and, from 


the lower curves, a new X. 


6. This process is continued until the solution 
converges. This will usually require three 
or four approximations. 
The solution can be somewhat shortened by replotting 
the data as in Figure J-2. Here flux time has been eliminated 
between the two curves in Figure J-1. In addition, the value 


of Nep/Noe plotted for each value of X has been modified by 
O 
the value of X so that Npp/Nosr can be read directly. In 


finding X by Figure J~2, Nep/No « is found for the reference 
design and establishes a vertical line which needs no fur- 
ther modification. Along this line, ICRpR can be used to 
find an initial estimate of X. This X is then used to modi-~ 
fy ICR as in step (5) above, and the process repeated until 
convergence. Figure J-1 is a somewhat clearer representa- 
tion of the logical steps in the determination of X, but 
Figure J-2 will give slightly quicker convergence than 
Figure J-1, eliminates the modification process on Nep/No 5 ; 


and concentrates the data in one set of curves. 
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